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CHAPTER  I 


INTRODUCTION 


General 


When  the  design  of  protective  structures  to  resist  blast  and  other  close-in  nuclear  weapons  effects  is  considered,  it  is 
only  logical  to  realize  that  simultaneously  large  amoums  of  initial  and  fallout  radiation  may  be  inevitable  realities.  The 
shielding  that  is  required  to  provide  adequate  resistance  to  close-in  radiological  effects,  and  the  structural  strength 
required  to  properly  resist  the  blast  effects  in  these  regions,  suggest  buried  structures  as  the  most  feasible  solution.  It  is 
observed  that  the  buried  structure,  in  addition  to  providing  resistance  to  these  effects,  possesses  a  hardness  which  is 
applicable  to  resisting  other  effects,  for  example  those  from  the  immediate  thermal  pulse  and  the  effects  of  later  fires. 

Laymen,  and  engineers  alike,  have  for  years  concluded  that  such  shelters  are  beyond  the  reach  of  economical 
programs.  This  in  fact  is  reflected  in  current  national  OCD  policy  and  the  National  Shelter  Program.  However,  such 
reasoning  is  usually  not  based  on  studies  which  consider  true  minimum  cost  structure*.  In  general,  the  structures  on  which 
these  conclusions  are  based  are  traditional  in  concept.  These  traditional  structures  in  general  grew  from  application  of 
traditional  design  procedures  to  situations  in  which  loadings  are  of  one  or  two  orders  of  magnitude  greater  than  traditional. 

Traditional  Configurations 

Traditional  building  configurations  or  shapes  were  originally  conceived  and  generated  to  hou^e  man  in  everyday 
fractions;  and  were  called  upon,  structurally,  to  resist  only  naturai  and  self  induced  loads  of  moderate  magnitude.  Most 
of  the  engineering  technology  developed  for  traditional  cr  conventional  design  was  directed  toward  the  development  of 
procedures  and  methodologies  to  adequately  assu.'e  the  designer  that  this  structure,  which  for  example  was  composed  of 
beams  and  slabs,  would  properly  resist  these  conventional  loads.  The  magnitude  of  load  rarely  affected  very  greatly  the 
choice  of  structural  type.  When  overpressures  from  nuclear  bursts  are  considered,  and  then  the  existing  design  procedures 
are  applied  to  design  traditional  type  structures  for  these  higher  loadings,  there  is  little  wonder  that  expensive  designs  result. 

Design  Economy 

Two  aspects  (among  others)  are  present  in  the  design  of  all  civil  engineering  structures.  They  are  the  selection  of 
configuration,  and  the  determination  of  size.  The  configuration,  or  shape,  is  a  qualitative  matter  v/hich  describes  in 
general  terms  the  nature  of  the  structure  such  as  a  flat  slab,  beam  slab,  concrete  joist,  or  dome  type.  The  other  aspect  is 
size  of  the  structural  components  which  are  contained  in  the  shape.  These  are  aspects  such  as  thicknesses,  amounts  and 
location  of  reinforcement,  etc.  In  the  traditional  civil  engineering  structure,  the  shape  decision  has  generally  been 
somewhat  arbitrary.  The  size  decision  , within  the  shape,  has  been  determined  by  more  or  less  rigorous  calculations. 

This  procedure  results  in  reasonably  economical  solutions  if  loads  and  spans  ore  small.  In  the  case  of  design  for  blast 
resistance,  the  loads  are  significantly  higher  than  any  considered  previously.  The  spans  may  be  shortened,  but  usually 
at  a  loss  of  functional  value.  Thus  it  is  not  surprising  to  discover  that  selection  of  configuration  can  no  longer  be  as 
arbitrary . 

In  summary,  if  the  structural  configuration  is  selected  without  regard  to  the  magnitude  and  type  of  loading,  the 
resulting  size  determinations  will  invariably  lead  to  the  requirement  of  a  great  amount  of  material,  and  therefore  higher 
cost.  If,  however,  structural  configurations  are  chosen  out  of  proper  respect  for  loads,  the  resulting  sizes  will  be 
reasonable  and  hence  material  quantities  which  are  related  to  cost  will  be  reduced  and  held  near  a  minimum. 

Efficient  Protective  Structures  Configurations 

The  most  efficient  forms  that  man  knows  to  resist  pressures  are  doubly-curved  shell  structures.  A  shell  structure 
achieves  its  efficiency  ly  resisting  loads  primarily  through  the  development  of  direct  stress.  Examples  of  these  types  are 
the  spherical  and  parabolic  domes  that  are  common  in  civil  engineering  solutions  to  the  problem  of  spanning  large  arenas 
with  a  minimum  of  material.  The  shell  of  an  egg  is  another  common  example  of  this  form  which  exhibits  the  same  type  of 
structural  efficiency,  Amonn  the  classes  of  doubly-curved  shells,  which  ore  available,  are  some  which  exhibit  truly  the 
ultimate  in  structural  efficiency.  These  are  the  constant-stress  funicular  shells.  The  spherical  pressure  vessel  is  an 
example  of  this  type. 

Constant  stress  funicular  shells  resist  pressure  loads  by  jither  developing  uniform  tensile  or  uniform  compressive 
stresses  depending  on  the  particular  configuration.  The  concrete  shell  works  best  in  compression;  the  steel  shell  in  tension. 
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Unfortunately,  compressive  funicular  shells  which  have  low  thickness  to  curvature  ratios  tend  to  buckle  on  overload.  As 
a  general  statement,  a  shell  structure  (or  any  structure)  in  a  state  of  compression  may  buckle  catastrophically  if  certain 
combinations  of  geometric  and  material  parameters  are  ill  chosen.  Because  the  calculation  or  prediction  of  this  tendency 
is  a  complex  procedure  for  all  but  the  simplest  of  structures;  and  at  best  is  an  approximation,  it  is  probably  better  to  avoid 
the  problem  entirely  by  insisting  that  tensile  funicular  shells  be  used.  These  structures  cannot  buckle. 

The  Contract 


This  report  is  the  result  of  a  contract,  a  quote  from  which  follows.  It  investigates  the  feasibil  ty  of  the  application 
of  these  funicular  structures  to  the  problem  of  sheltering  people  from  dynamic  overpressures  caused  by  nearby  nuclear  or 
conventional  bursts. 

"A.  The  contractor,  in  consultation  and  cooperation  with  the  Government,  shall  furnish  ali 
engineering,  labor,  equipment,  tools,  materials,  supplies,  facilities,  and  services  necessary 
for  a  feasibility  study  relating  to  optimizing  shelter  design.  The  work  and  services  shall 
pertain  to  the  analysis  and  design  of  flexible  yielding  membrane  elements  of  a  shelter  to  resist 
normal  dynamic  effects  not  unlike  those  which  may  result  at  the  soil -structure  interface  as  a 
result  of  a  nuclear  blast. 

B.  The  general  areas  of  investigation  shall  include,  but  not  be  limited  to  the  following: 

1 .  Investigate  the  theoretical  prediction  of  the  configuration  of  a  yielding  membrane 
and  determine  its  application  to  the  shelter. 

2.  Perform  certain  loading  simulator  studies  to  corroborate  the  intuitive  foct  that  yielding 
buried  structures  are  efficient  structural  systems. 

3.  Extend  the  theory  of  studies  involving  the  investigations  of  the  membrane  supported  on 
yielding  boundaries. 

4.  Determine  the  feasibility  of  future  possible  exploitation  in  this  area," 


2 


CHAPTER  il 


General 


SHELTER  APPLICATIONS  AND  CONSTRUCTION  TECHNIQUES 


In  the  course  of  introducing  new  concepts,  such  as  that  of  the  yielding  membrane,  it  is  often  desirable  (from  a 
presentation  standpoint)  to  indicate  the  ultimate  uses  of  the  product  before  the  involved  supporting  calculations  and  data 
are  presented.  This  pattern  is  used  here;  therefore,  the  material  in  this  Chapter  serves  in  part  os  illustration  of  an 
application  of  the  concept.  Models  are  used  as  examples  of  possible  shelter  applications. 

The  f:rst  example  illustrates  the  rather  limited  application  of  a  discontinuous  circular  dish  membrane  as  the  roof 
structure  of  a  community  shelter  which  is  designed  to  resist  55  psi  blast  overpressure  from  a  5  MT  nuclear  weapon  which  ;< 
detonated  as  a  surface  burst.  The  second  is  on  attempt  at  approaching  the  ultimate  in  efficiency  of  this  application. 

This  approach  involves  a  multistory  cubicle  structure  in  which  the  yielding  membrane  is  continuous  and  completely 
encloses  the  cubicle  volume. 

Circular  Membrane  Structure 


The  model  shown  in  Figures  1  through  6  illustrates  the  application  of  the  concept  through  the  use  of  a  dished 
membrane  as  a  roof  structure.  Although  blast  simulator  studies  (Chapter  III)  indicate  that  the  preliminary  dishing  is  in 
fact  unnecessary,  from  a  structural  standpoint,  the  application  .hown  is  a  conservative  one  and  is  to  be  recommended 
until  more  test  and  analytical  data  are  obtained. 

The  steel  cap  or  membrane  carries  blast  overpressure  loadings  primarily  by  the  development  of  direct  uniform 
tensile  stress.  Little  or  no  bending  exists;  hence  little  or  no  shear  is  present.  In  the  application  shown,  the  membrane 
is  not  continuous.  That  is,  it  is  a  series  of  individually  dished  elements.  This  break  in  continuity  requires  the  use  of 
boundary  arch  structures  which  are  the  reinforced-concrete  edge  rings.  They  are,  in  fact,  concrete  arches  which  lie  in 
a  horizontal  plane.  They  eventually  transmit  the  vertical  component  of  the  load  to  the  supporting  wall  structure. 

Note  the  visual  expression  of  structural  efficiency  which  Is  present  when  steel  is  used  in  tension  and  concrete  in  compression. 

To  provide  added  toughness  in  design  concept,  the  modular  system  is  used.  Nature  has  used  the  same  system,  that 
of  duplication  of  self  sustaining  elements,  in  many  naturally  occurring  designs;  in  order  to  assure  survival  of  the  function 
of  the  overall  product  against  attacks  by  facets  of  hostile  environments ,  These  modular  units  are  30  feet  in  diameter  and 
conservatively  use  3/8  inch  steel  membrane  roofs.  Each  unit  is  an  independent  structure  which  does  not  rely  upon  the 
strength  of  adjacent  units  to  provide  a  reaction  or  to  contribute  to  the  support  of  external  forces.  The  particular  advantage 
of  this  indeoendent  action  is  that  the  other  elements  of  structure  would  not  collapse  if  one  unit  failed.  It  must  be  noted 
here  that  large  structures  would  be  subjected  to  unsymmetrical  loading  caused  by  the  transient  nature  of  the  blast  wave, 
which  would  tend  to  shear  each  unit  from  the  other  vertically.  The  design  suggested  adjusts  to  this  shearing  tendency. 

The  prototype,  which  the  model  illustrates,  would  lend  itself  readily  to  precast  and  prefabricated  construe  ion 
methods.  The  entire  wall  structure  con  be  constructed  from  one  basic  snape.  Precast  concrete  piers  may  be  sat  in  ploce. 

Roof  cups  may  be  field  welded  to  attachment  points.  A  variation  of  this  concept  includes  a  membrane  floor  of  the  same 
configuration,  and  material,  as  the  roof.  An  intern, oi  floor  system  would  be  required  as  well  as  other  changes  in 
construction  technique.  Cost  estimates  for  this  structure  were  reported  previously  on  other  contracts.* 

Rectangulor  Membrane  Structure 

if  opt!mization  of  efficiency  of  configuration  and  application  to  the  total  shelter  is  attempted  In  concept,,  one 
solution  which  results  is  that  presented  in  Figures  7  through  ?.  These  a. •  views  of  a  model  which  it  nearly  the  shape  of  a 
cube.  The  cube  is  the  rectilinear  solid  which  has  the  maximum  interior  volume  for  tho  minimum  exterior  surface.  Secusne 
blost  overpressures  in  shelters  are  hopefully  resisted  oniy  a?  the  exterior  suit  ores,  the  structure  that  incorporates  o  minimum 
of  these  surfaces  uses  a  minimum  of  material  within  the  configuration  or  shape  cluirificotion .  Such  is  o  fundamental 
requirement  of  material  cost  minimization. 


*  Final  Report,  "Local  Civil  Defense  System*,"  Contract  OCD-OS-62-232,  University  of  Arizona,  June,  1964 
and  Final  Report,  "Cost  Studies  in  Protective  Construction  Systems,"  Subcontract  Institute  for  Defense  Analy*-n, 
1 38—4,  Jonuory,  1965. 
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Figure  3  .  Crow-reef  lor  of  -community  shelter  module 
showing  disn-type  membrane  roof  structure. 


Figure  4  .  C  lore -up  of  construction  procedure 


Figure  5  .  Construe  Hon  procedure  onci 
sloping  entrance 


Figjre  6  •  Precast  column  .■.hicii  assembi*s 
to  form  c!l  inferior  *ol!s 
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Figure  7  (Right) .  A  cutaway  view 
of  a  model  which  illustrates  the 
concept  of  the  use  of  rectangular 
yielding  membranes  as  exterior 
structural  parts  of  a  shelter  complex. 
This  model  uses  a  minimum  of 
exterior  skin  to  a  maximum  of 
interior  volume.  The  "egg  crate" 
walks  and  floors  provide  the  needed 
rigidity  to  keep  the  membrane  in 
place. 


Figure  8  (Left).  The  membranes  take  a  dishea  Hrm  c 
indicated  in  this  photograph  when  they  are  subjected 
to  external  pressure.  This  dishing  process  promotes 
favorable  soil -structure  interaction  behavior. 


Figure  9  (Above).  Requirements  of  flexibility  in 
connections  are  crudely  illustrated  by  this  detail. 

It  is  necessary  ta  use  flexible  links  between  the  tva;ic 
shelter  and  its  somewhat  disconnected  entro'sce. 
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To  further  optimize  the  use  of  material,  the  exterior  skin  of  this  shelter  is  a  continuous  yielding  membrane.  As  in 
the  previous  example,  dished  elements  are  used  even  though  tests  indicate  such  preliminary  dishing  is  unnecessary.  The 
advantages  of  continuity  are  apparent  in  the  savings  produced  by  the  lock  of  reinforced  concrete  supporting  edge  rings. 
The  interior  walls  and  floors,  which  ore  of  reinforced  concrete,  transfer  the  thrust  through  the  entire  shelter  where  it  is 
equilibrated  by  the  loads  on  the  opposite  sides.  This  shelter  overcomes,  in  port,  many  of  the  structural  disadvantages  of 
large  one-story  structures.  Aside  from  using  a  minimum  of  exterior  surface,  it  resists  the  "slicing  up"  shears  that  occur  in 
most  large  structures  as  the  blast  wave  travels  across  the  shelter.  It  is  noted  in  passing  that  these  can  be  extremely  large 
but  they  are  not  usually  considered  in  normal  routine  design  practice. 

The  structural  shell  as  illustrated  possesses  reserve  strength  to  resist  all  conceivable  ground  motions  associated  with 
any  overpressures  for  which  it  is  designed.  Interior  details,  of  course,  must  be  designed  in  such  a  way  that  the  creation 
of  interior  missiles  and  associated  damage  to  occupants  is  prevented  * 

Other  Applications 

No  doubt,  in  the  review  of  the  applications  presented  in  this  Chapter,  other  applications  of  the  membrane 
cor.cept  are  indicated.  These,  in  fact,  do  exist  and  include  such  items  as  blast  doors,  end  walls  to  tube  structures  and 
above-ground  protective  enclosures  for  supplies;  and  other  items  which  are  unaffected  by  such  factors  as  electromagnetic 
pulse,  initial  nuciear  radiation,  and  fallout  gamma  radiation.  Extremely  thin  membranes  can  resist  very  high  overpressure 
in  such  applications  but  they  must  be  shielded  from  missile  penetration,  in  most  cases,  it  is  not  necessary  for  such 
membranes  to  be  pre-formed.  Their  ultimate  resistance  to  blast  overpressures  is  unaltered  by  leaving  them  initially  flat. 
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CHAPTER  Ml 


ENGINEERING  ANALYSES 


Introduction 

Intuitively,  it  appears  thot  yielding  membranes  are  technically  sound  protective  structure  components.  To  support 
this  intuitive  reasoning,  it  is  appropriate  to  consider  certain  supporting  calculations  and  analyses  of  an  engineering 
nature.  In  addition  to  the  analysis  of  the  continuum,  attention  must  be  given  to  the  boundaries  where  the  membranes 
terminate. 

With  respect  to  these  problems  the  details  of  the  historical  development  of  analysis  of  these  structures,  together 
with  specific  developments  for  the  instant  shapes,  are  presented  at  length  in  Appendices  A  and  D  to  this  report.  A 
summary  of  some  of  the  more  essentio!  engineering  features,  which  are  in  port  abstracted  from  this  analysis,  are  included 
in  this  Chapter. 

General  Membrane  Theory 

Flexible  buried  structures.-  because  of  yielding  characteristics  which  produce  negative  settlement  ratios,  offer  the 
ultimate  in  economy  in  view  of  the  way  in  which  they  force  the  soil  to  resist  the  overload.  The  most  efficient  flexible 
structure  is  that  which  simultaneously  yields  under  constant  stress  at  every  point  in  its  plane.  The  behavior  of  such  o 
structure  may  be  predicted  in  advance  by  an  inverse  solution  of  the  differential  equations  for  stress  in  shell  structures 
under  normal  pressure  loadings.  These  structures  by  definition  ore  called  funicular. 

For  an  example  of  the 
structural  efficiency  of  a  system 
such  as  this,  consider  the  following 
simple  comparisons.  These  structures 
ere  intentionally  not  buried  for 
simplicity  of  presentation. 

Membrane  Analysis  of  a  Thin  Plate 

For  a  comparison  consider  a 
previously  designed  one-way 
fiat  slob. 

Where  under  a  standard 
design  with  pQ  =  50  psi  the 
resulting  section  is  as  follows: 

Now  consider  the  some 
span  covered  with  a  thin  steei 
plate. 

Assume  that  when  the  load 
is  applied  the  plate  will  yield  into 
a  circular  arc.  Consider  a  free 
bcd>  of  the  loaded  section,  as 
shown  on  the  next  page. 
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By  static  Equilibrium; 

£f  =0 

V 

PC(L)  "  2(10,000  )  =  0 

„  _  10,000  „  „  _  10,000 

Po - or  R  '  P 

o 

For  this  particular  case  of  loading: 

Pq  =  50  psi  .  .  R  = — ^ — =  zOO  in. 

This  determines  the  Radius  as  a  function  of  the  load  only, 
independent  of  the  length  (L).  This  condition  in  itself  is 
insufficient  since  no  consideration  is  given  to  the  percent¬ 
age  elongation.  To  determine  this  percentage,  consider 
the  equation  of  the  triangle  bounded  by  R  and  R-h. 

2  2  L  ^ 

(R-h)  -  R  -  (j) 


!>=■*-  Vr2-<5>2 

h  for  this  cose  is  given  by  h  =  200  -  [  (200)^ 
Calculation  of  the  percentage  of  elongation: 


At  yield  T  =  f  t  =  40,000  (1/4)  -  10,000  Ibs./iri. 


200-  178  =  22" 


%  e  = 


G-L 

L 


where  G  is  the  arc  length 
of  the  membrane 


-1  90 

Q  =  R6  =  (200)  2  tan  (-^  =  187.2' 


%  e  = 


187.2  -  180  _  7.2 
180  "’80 


=  0.04  =  4%  .\  C.K . 


Thus  it  is  seen  tiiat  a  1/4"  steel  ijjembrane  is  capable  of  resisting  as  great  a  load  cs  a  19"  reinforced  concrete  slab  which 
is  reinforced  at  a  rate  of  4.08  in  /fr.  The  plate  contains  3.0  in^/ff*  —  less  steel  than  in  the  reinforced  slob.  The  strain 
of  4%  is  less  than  the  ultimate  uniaxial  strain  capacity  of  most  structural  steel  plate. 

For  increased  efficiency  of  steel  membrane,  it  may  be  used  in  a  biaxial  state  of  stress  —  such  os  that  found  in  a 
circular  diaphragm,  cn  example  of  which  follows. 


Two-Way  Circular  Membrane 

Consider  a  circular  plate  of  diameter  L 
subjected  to  a  load  of  p^  psi  and  clamped 
around  the  circumference,  as  shown  at  the 
right. 
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From  static  equilibrium: 

F  -0 
v 

2  2 

PQtr  =  10,000  (^)Tr^) 

.  D  20,000  .  20,000 
••  po=“i —  or  R"— f“ 

o 

For  P  =  50  psi,  R  =  400  in , 
o  i 

The  same  examples  may  be  solved 
a  more  general  approach  This  a 

the  application  of  the  general  membrane  theory  of  shell 
structures  to  the  situation  presented.  An  introduction  of 
this  approach  to  design  was  made  at  the  Symposium  on  Shell 
Research,  Deift,  The  Netherlands,  August  30,  1961,*  The 
application  at  that  time  was  directed  toward  the  "Configuration  of 
Shell  Structures  for  Optimum  Stress."  Basically,  the  approach  involves 

the  initial  assignment  of  a  given  final  srress  stete,  such  as  that  of  constant  stress.  The  search  is  then  made  for  the  shell 
structure  v/hich  exhibits  this  find  state  of  stress  under  a  previously  assigned  normal  pressure  loading. 

Consider  a  free-body  of  an  element  of  a  shell,  as  shown  above.  From  the  equilibrium  conditions,  if 
-  N2  *  +  S  (constant),  and  ^2  =  ^21  ~  ^ 


fhen  Rj  +  J<2  "  S 

Note  that  in  this  equation,  Rj  and  R2  are  the  principal  radii  of  curvature  of  the  final  deflected  surface.  P  is  the 
normal  pressure  on  the  surface  end  S  is  the  membrane  tension  in  dimensions  of  force  per  unit  length. 

An  illustration  of  the  application  of  this  equation  is  now  made  with  reference  to  the  previous  two  examples. 


*  H.  P.  Harrenstien,  "Configuration  of  Shell  Structures  for  Optimum  Stresses,"  Proceedings  of  the  Symposium  on 
Shell  Research,  Delft,  The  Netherlands,  1961, 
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Second  Example: 


Obviously,  if  is  as  single  to  apply  the  free-body  diagram  approach  as  it  is  to  apply  the  general  theory,  however 
for  mere  complex  problems,  the  general  theory  must  be  used  in  conjunction  with  numerical  solutions  on  a  digital 
computer.  Appendix  A  presents  this  approach,  the  results  of  which  are  abstracted  in  the  following  material. 

General  Theory  of  Funicular  Shells 

A  shell  structure  may  be  defined  as  a  materialization  of  a  curved  surface  in  space.  In  general,  the  structure  of  a 
structural  component  so  formed  carries  its  loads  primarily  by  direct  stress.  By  this  process,  such  structures  make  maximum 
use  of  the  material  from  which  they  are  formed  in  resisting  applied  loads  (see  Appendix  B). 

As  stated  before,  among  the  classes  of  sheli  structures  that  are  available  for  shelter  application,  there  exists  those 
which  initially  exhibit  uniform  direct  stress  characteristics  under  certain  specified  loadings.  These  types  of  shell 
structures  possess  the  maximum  possible  structural  efficiency  that  a  two-dimensioial  structure  is  capable  of  providing. 

A  structure  which  is  of  one  sheet  and  which  exhibits  uniform  stress  characteristic;  under  normal  pressure  loading  is 
defined  as  a  funicular  shell.  The  particular  type  of  such  shell  structures  that  are  formed  when  a  thin  steel  plate  yields 
with  "constant"  stress  under  application  of  a  distributed  normal  pressure  is  one  such  shape  and  is  the  subject  at  hand. 

The  pressure  need  not  be  uniform  for  the  funicular  concept  to  be  present. 

The  exact  equation  of  curvature  derived  from  differential  geometry  for  a  function  z  =  f(x,y)  is: 
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Solutions  and  Design  Curves 

This  expression  may  be  written  in  finite  difference 
form  and  solved  by  iteration  techniques  on  a  digital 
computer  for  certain  membrane  shapes  and  specific  edge 
conditions,,  The  types  considered  are  shown  in  Figure  10. 
The  non-dimensionalized  PD/S  versus  Zg/D  curves,  which 
are  shown  in  Figure  1  l,are  results  from  the  computer  study. 
Only  the  curves  showing  the  center  deflection  are  given, 
but  these  are  the  most  important  as  far  as  design  is  con¬ 
cerned.  These  curves  are  based  on  the  behavior  of  a 
rigid-plastic  material .  However,  to  use  these  curves 
with  any  other  type  of  material  the  only  additional 
information  required  is  the  appropriate  stress-strain  curve 
for  the  material . 
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The  PD/S  vs.  Zq/D  curves  are  based  on  the  average 
stress  and  average  strain  across  the  center  of  the  membrane. 
It  is  known  that  the  strains  are  not  uniform  over  a 
deflected  membrane  surface  (2,3,4).  However,  as  far  as 
vertical  deflections  are  concerned,  the  assumptions  of 
unifo  m  stresses  and  strains  appear  justified. 

The  only  regions  where  this  assumption  leads  to 
appreciable  errors  is  in  the  corners  of  rectangular 
membranes.  If  reasonable  care  is  taken  during  the  con¬ 
struction  to  insure  proper  full-strength  weids  and  if  the 
design  strain  is  reasonable  (less  than  2.5%)  the  yielding 
membrane  structured  element  should  serve  quite  well. 

Either  the  circular  or  rectangular  problems  could  be 
programmed  for  the  computer  with  a  non-uniform  lateral 
pressure.  Soon,  it  may  be  possible  to  predict  the 
attenuation  of  pressure  on  a  yielding 
structural  element  and  the  resulting 
pressure  distribution.  However, 
in  working  with  yielding  elements, 
they  can  be  designed  as  if  they 
were  to  be  subjected  to  the  full 
uniform  lateral  pressure.  The 
yielding  characteristics  force  the 
surface  to  take  the  shape  it  must 
assume. 
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Clomped  Edge 


Figure 


Yielding  Beam 

10.  Plan  Views  of  Shapes  Considered 
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Although  certain  metals, 
especially  mild  steel,  have  very 
large  plastic  elongation  propenies 
on  uniaxial  tests— -sometimes 
greater  than  30%  strain— this  doe: 
not  mean  the  material  will  admit 
such  large  strains  under  biaxial 
conditions.  In  fact,  most  of  the 
common  yielding  materials  will 
not  admit  average  strains  greater 
than  9  to  10  percent  even  in  a 
membrane  state  of  stress.  Since  a 
true  membrane  state  of  stress  may  be 
difficult,  if  not  impossible  to 
realize  in  actual  construction,  a 


PD 

s 


/  o 


Figure  1 1 .  PD/S  versus  z  J D  Curves 
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maximum  design  strain  of  2.5%  is  recommended.  This  strain  corresponds  to  a  z^D  ratio  of  about  0.10. 


To  illustrate  the  use  of  these  design  curves,  consider  the  circular  membrane  which  was  solved  previously.  In  that 
example 


For  this  case 


r  -  v  =  ou  psi 
so  r 

i_ 

-  - ^ 

*v* , 

1  ••  • 

D  =  15  ft  =  180  in. 

»•  •  ■* 

S  =  10,000  Ib/in. 

;  .1; 

0  *• 

PD_(50)(180) 

V  -  * 

"Shock  Isolated"  Floor 

l »’ 

S  10,000  w 

:'v 

j  f  1  — t 

0\ 

•  • ; 

— 7“ 

: 

W  •  C>  - -  - - - 1 - -  “  W. 

1%.  The  center  deflection,  z£,  is  then  0.057  D,  or  10.3  in. 

The  ultimate  strength  of  this  membrane  may  be  easily  determined  by 
entering  Figure  11  with  a  maximum  biaxial  strain  requirement  of  2.5%. 
For  this  strain  PD/S  =  1 .5  and  z JD  =  0.097.  If  PD/S  =  1 .5,  then 


<— -  Yielding  Membrane 
Figure  12.  Yielding  Membrane  Shelter 


1.5  S  _  1.5  (10,000)  _ 

p  d - rsc - 83ps* 


The  center  deflection,  zc,  for  this  condition  is  zc  -  0.097  (180)  -  17.5  in.  This  example  demonstrates  the  remarkable 
reserve  strength  of  these  elements. 


To  achieve  the  greatest  economy  and  overall  toughness  of  the  shelter,  it  is  suggested  that  the  same  type  of  membrane 
be  used  on  fhe  floor  as  on  the  roof.  Figure  12  above,  illustrates  this  concept. 


Usually  the  center  deflection  to  span  ratio  will  be  the  governing  design  factor  but,  also  a  check  should  be  made  to 
insure  against  an  excessive  pressure  increase  in  the  structure  which  may  be  induced  by  the  sudden  deflection  on  the  roof. 


This  "back  pressure" 
should  not  be  greater  than 
4  to  5  psi.  The  Lovelace 
Foundation  indicates  that 
this  is  the  threshold  of  the 
eardrum  damage  region. 

The  back  pressure  curve 
which  is  shown  in  Figure  13 
is  for  a  circular  structure 
but  will  work  well  for  | 

square  areas.  If  used  for  £ 
other  rectangular  shapes,  £ 
the  actual  pressure  would  “ 
be  greater  than  the  value 
from  the  graph  resulting  in 
non-conservative  answers. 

The  whole  problem  of 
back  pressure  can  be 
ignored  if  the  membrane 
has  an  initial  "dish." 
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Membranes  Supported  by  Yielding  Beams 
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If  yielding  membranes  were  used  in  fhe 
design  of  blast  shelters,  it  might  be  advantageous 
to  use  yielding  beams  as  ribs  across  the  membrane 
to  decrease  the  maximum  deflections.  As  with 
the  yielding  membrane  itself,  the  force  in  fhe 
yielding  beam  would  have  to  be  resisted  in  some 
manner.  In  the  following  figures  the  behavior 
of  such  reinforced  membranes  is  indicated. 


The  PD/S  versus  zc/D  curves  (Figures  14, 

15,  and  16)  are  non-dimensionclized  pressure 
versus  deflection  curves  for  the  points  of  maximum 
deflection  of  the  membrane  and  the  center  points 
on  the  beams  for  the  conditions  of  edge  constraint. 
The  numbers  inside  the  circles  are  values  of  fhe 
ratio  F/SD  where: 

F  =  strength  of  beam  and  membrane 
S  =  membrane  strength 
D  =  short  span  distance 

The  subscripts  refer  to  the  locations  of  the  point. 
For  example,  (T)  ^  ^  refers  to  the  PD/S  vs.  z/D 
curve  for  the  center* point  of  the  membrane,  when 
the  ratio  of  beam  strength  to  the  product  of  the 
membrane  strength  and  the  span  is  one.  These 
graphs  are  for  symmetrical  cases,  i.e.,  it  is 
assumed  that  the  conditions  on  both  sides  of  the 
supporting  beams  are  the  same  (see  Figures  18  and 

19). 

Boundary  Conditions 

The  yielding  membrane,  to  be  effective, 
must  be  supported  by  bounding  structural 
elements  which  are  capable  of  absorbing 
the  full  thrust  of  the  material,  at  yield. 

For  the  type  of  situation  shown  in  the  first 
example  of  Chapter  II,  the  concrete  arches 
are  well  suited  to  provide  the  necessary 
support.  For  the  second  example  shown  in 
Chapter  II,  other  details  must  be  considered. 

The  following  considers,  in  detail,  examples 
of  preliminary  design  solutions  to  the  enumerated 
types  of  boundary  support  problems  (Figure  17): 

1 .  Yielding  membranes  supported  by 
concrete  arch  rings 

2.  Yielding  membranes  supported  by 
straight  edge  beams 

3.  Yielding  membranes  supported  by 
curved  edge  beams 

4.  Yielding  membranes  supported  by 
in  plane  beams 

5.  Yielding  membranes  supported  by 
membrane  elements 
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Yielding  Membranes  Supported  bv  Concrete  Arch  Rings  (Figure  17a).  For  on  example  of  a  supporting  reinforced 
concrete  arch  ring  for  the  circular  memorane,  see  the  design  on  page  19.  Assume  that  the  span,  L,  is  30  feet,  a  con¬ 
siderable  span  for  a  1/4  inch  roof  structure  which  supports  50  psi.  The  edge  ring  requirement  may  be  determined 
approximately  as  follows: 


Consider  the  section  loaded  as  shown  by  the  full 
yield  strength  of  the  membrane.  Note  that  at 
initial  yield,  the  membrane  force  is  horizontal. 

\’o:  preliminary  calculation  purposes,  assume  a 
square  arch  of  cross  section  b  x  b  and  reinforced  by 
mi  Id  steel  at  a  percentage  of  4%.  A  free  body  diagran 


E  Fx  =  0  yield: 

2P  =  (10,000)  (30)  (12) 
P  =  1.8  x  106  lb. 


It  is  assumed  that  .he  ultimate  concrete  strength  is  3,000  psi  and  that  the  yield  strength  of  the  reinforcing  steel  is 
40,000  psi.  Equation  (19-7)  of  the  1963  code  of  the  American  Concrete  Institute  applies  to  this  situation  —  that  is,  a 
'olumn  which  is  under  direct  compression. 


✓here 


P  -  0.85  f1  (A  -  A  )  +  A  f 
0  e  g  St  st  y 


is  the  ultimate  axial  load 

is  the  ultimate  concrete  strength 

is  the  gross  section 

is  the  area  of  the  steel,  and 

is  the  yield  strength  of  the  steel. 


Substituting 


A  =  0.04  b 
st 


and 


A  =b2, 
9 


we  heve 


P  =  0.85  f1  (b2  -  0.04  b2)  +  0.40  b2  f 
o  c  y 

Pq  =  0,799  P  b2  +  0,04  f  b2  =  b2  [.799  (3,000)  +  .04  (40,000)] 
=  3997  b2 


Equating  to  P: 

1.8  x  106  =  3997  b2 
b2  =  450.34 


b  »  21.7  in. 

A  =  0.04  (b2)  =  0.04  (450.34)  =  18  in.2 
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Use,  for  preliminary  purposes,  the  following  section: 


5/16"  continuous  fillet  welds 


Me 


2000  (3.25)  &25)  6 


40,000  = 


21125 


t 


t  =  0.53 
f  =  0.73  in. 


5,000  Ib/in. 


Use  4  x  4  x  3/4"  angles  with  t  =  0.75  in. 

Determination  of  area  of  reinforcing  steel: 

2 

Two  4  x  4  x  3/4  angles  have  an  area  of  10.88  in  . 

.  ...  2  .2 

The  required  is  18  in  ,  therefore  7.12  in  mus*  be  added  as  reinforcement. 

2 

If  12  bars  are  used,  per  bar  is  0.595  in  . 

2 

Use  12  *7  V  0.60  in  per  bar. 


5/16"  continuous  fillet  weld 
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Yielding  Membrane  Supported  by  Straight  Edge  Beams  (Figure  17b).  It  is  logical  to  consider  the  method  of 
supporting  rectangular  membranes  at  their  outer  edges  by  initially  straight  steel  beams  which,  on  loading,  yield  inward 
much  in  the  same  manner  as  the  yielding  stiffening  ribs,discussed  previously.  If  the  design  of  such  are  considered,  it 
becomes  obvious  that  the  amount  of  material  required  is  excessive  and  such  o  support  <s  impractical.  This  Is  demonstrated 
in  the  following  development. 


pninnEE 


rfr$ 

rmr 


S  Ib/in. 


Consider  a  straight  edge  member  which  is  loaded  by 
a  uniform  membrane  tension.  Let  the  span  be  L  between 
supports  and  the  tension  be  S  in  ib.  per  in.  The  situation 
is  shown  at  the  right.  Under  action  of  force  S,  the  beam 
deflects  in  a  circular  arc  of  radius,  R,  to  a  maximum 
center  deflection,  h.  The  beam  in  this  configuration  carries  a  maximum  axial  load  of  P  Ib.  A  free-body  diagram  of  the 
deflected  member  is  shown  as  follows: 


R 


P  =  SR 


LFy  =  0  yields 

2P(^)  =  $L 


For  a  given  steel  beam  with  an  area.  A,  and  a  yield  stress,  f^, 

P  =  f  A 

y 

on  substitution  and  simplification 


y 


Thus,  for  a  given  situation,  the  area  of  edge  beam  that  is  required  is  only  a  function  of  the  membrane  tension, 
the  yield  stress,  and  the  radius  of  curvature  of  the  deflected  shape.  Now,  S  is  the  membrane  yield  tension  which  is 
f  t  where  t  is  the  thickness  of  the  membrane.  Cn  substitution 

y 

f  I  R 

A  =  —  *  t  R 

y 


If  limits  are  established  for  the  axial  strain  in  the 
edge  beam,  say  about  10%,  then  R  has  a  limit  based  on 
this  strain.  To  determine  the  upper  limit  of  R,  which  is 
actually  the  minimum  value  of  R,  consider  the  geometry 
sho^n  at  the  right.  The  average  %  strain  is 


R  9  -  1/2 

C75 


(100) 


Now 

(100)  =  10 


Re- 1/2 


r  e  - 


11  L 
~W 
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t  =  1/4"  t  =  3/8' 


but 


R  arc  sin 


y±  ~ 

R  - 
L  11  L 


9  =  arc  sin  =  arc  sin  ^ 


2R 


20 


.  L  _  11  L  11  .L  , 
arCS,n2R=B  R=TO  2R 


arc  sin  :  - 


II 

115 


The  solution  to  this  transcendental 
equation  is  approximately  x  =  0.685.  If 
L/2R  =  0.685  for  10%  strain,  then 

R  =  2T^s5)  =  0,73  L 

and 

A  =  t  R  =  .73  t  L 

Solutions  to  this  equction  are  plotted 
on  the  figure  at  the  right. 


It  is  seen  that  even  vor  small  spans  of  around  10  ft,#  a  1/4  in.  membrane  requires  a  steel  edge  beam  of  21 .5  in. 
in  cross  section.  This  assumes,  of  course,  that  the  beam  is  initially  straight  and  that  it  yields  to  a  limited  maximum 
uniaxial  strain  of  10%. 


Yielding  Membranes  Supported  by  Curved  Edge  Beams  (Figure  17c).  If  the  edge  beams  are  ini'ially  curved  to 
radius,  R,  they  will  act  more  efficiently  when  tfie  membrane  yields;  because  they  are  not  limited  to  small  curvatures, 
by  strain  requirements.  For  this  case,  A  =  t  R  as  before.  This  equation,  for  various  membrane  thicknesses,  may  be 
plotted  as  shown  below. 


|t  is  observed  that  a  curved 
member  of  5  ft.  radius  and  area 
15  in^  would  provide  support  for  a 
1/4  in,  membrane.  This  produces 
a  slight  saving  in  material  over  the 
last  case,  but  is  much  more  difficult 
to  construct.  It  is  therefore  concluded 
tho*  such  support  is  impractical. 
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Yielding  Membranes  Supported  L>y  Reinforced  Concrete  Sheor  S^bs  (figure  1 7d) •  One  of  the  most  practical 
solutions  is  use  of  reinforced  concrete  slobs  in  cN  outside  panels  ond  cl  low  fheso  slabs  to  resist  the  membrane  thrust  of 
interior  elements  through  beam  action  as  deep  beams.  In  this  cose,  membrane  economies  will  on!/  be  cchieved  if  large 
numbers  of  repetitive  bays  are  used.  The  design  of  these  reinforced  concrete  slobs  (deop  beams)  is  fairly  routine  ar.d 
will  not  be  presented  here. 


Yielding  Membranes  Supported  by  Membrane  Elements  (Figur«  i  '<•).  A  method  of  support  that  seems  most 
likely  is  use  of  yielding  membranes  supported  by  membrane 


elemenrs.  At  the  corners  of  cubical  structures,  in  which 
yield  membranes  are  osea  as  the  exterior  structural  skin, 
membrane  elements  theinselves  may  be  used  to  "earner" 
the  thrust. 

For  this  case 

F  =  0  t-  f  cos  9  =  t .  f 

x  2  y  1  y 

t^  cos  0  =  fj 


if 

0  =  45° 

*2  "7737  =  MW  *1 

A  plot  of  the  relative  thicknesses  is  shown  above. 
For  1/4  in.  side  membranes,  0.35  in.  thick  comer 
membranes  are  required. 
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Membrane  Splices.  Practical  construction  considerations,  based  on  limited  sizes  of  available  sheet  steel, 
require  the  use  of  splices  on  membrane  plates.  Two  types  of  full  strength  splices  are  recommended.  They  are  shown 
as  follows: 


Simple  lap  splice 


Simple  butt  splice 


It  is  felt  that  either  of  these  splices  can  be  properly  designed  and  constructed  to  develop  the  full  yield  strength 
of  the  membrane  plates. 


Dynamic  tests  were  conducted  in  the  University  of  Arizona  blcet  simulator.  This  simulator  is  a  plane-wave 
generator  powered  by  a  hydrogen-oxygen  explosion.  The  blast  chamber,  shown  in  Figure  20,  is  a  300  gallon, 

8  ft*  by  2.5  ft  diameter  tank  mounted  vertically  on  rubber  bushings  to  a  heavy  concrete  base.  Access  to  the  chamber 
is  achieved  by  unbolting  the  top  section  of  the  tank  ond  swinging  the  bottom  section  and  the  base  around  a  pivot* 

The  2.5  ft.  x  2.5  ft.  soli  bin  is  then  exposed  for  the  placement  of  model  structures,  gouges,  and  sand.  There  is  also 
a  14  in.  diameter  access  Itolc  in  the  botTom  of  the  soil  bin  and  two  4  in.  diameter  access  holes  and  windows  in  the  too 
section. 

The  b^ttt  weve  is  caused  by  a  hydrogen-oxygen  explosion  detonated  by  an  electric  spark.  Predetermined 
quantities  of  hydrogen,  oxygen,  and  air  ore  meesured  in  rhe  three  auxiliary  tanks  on  the  side  of  the  chamber.  The 
oir  controls  the  rise  time  of  the  blast  wuve.  The  gases  are  fed  into  an  evacuated  plastic  bog  ar  the  top  of  the  tank. 

The  explosion  is  detonated  by  an  electric  spark  which  is  triggereJ  by  rhe  same  switch  that  sterts  the  recording  instruments. 
The  decay  time  is  controlled  by  adjustable  exhaust  valves  and  spacing  washers  between  the  chamber  sections.  The 
decoy  curve  is  exponential.  The  detonating  spark  and  the  gas  bag  are  centered  in  the  tank  to  minimize  dynamic 
imbalance  during  tests. 


The  overpressure  range  is  from  0  to  50  psi  with  variable  rise  times  from  less  than  one  millisecond  to  over  one- 
tenth  of  a  second,  and  decay  times  from  one-tenth  of  a  socond  ond  up.  The  blast  waves  could  be  controlled  to  within 
10  percent  from  rest  to  test.  The  instrumentation  and  recording  device*  included  two  Statham  pressure  transducers  (0  to 
50  p*i),  two  Tectronix  dual  beam,  oscillators  with  cameras,  one  six-channel  Brush  recorder,  and  one  two-channel 


Sanborn  recorder. 

The  blast  simulator 
has  been  used  to  test 
yielding  membrane  models 
above  and  below  ground. 
The  test  models  were  3.5 
to  4.0  inch  cyiinders  which 
were  constructed  so  that 
the  edges  held  c  membrane 
tightly  clamped  (Figure  21). 
The  investigations  were 
made  to  determine  the 
effects  of  depths  of  burial 
ond  structural  flexibility 


on  the  percentage  of  load 
carried  by  the  yielding 
elements. 

From  the  tests  made 
in  the  blast  chamber,  some 
insight  has  beer,  gained 
into  the  amount  of 
aitenuatioi:  of  overpressure 
which  is  caused  by  soil 
cover.  These  tests  showed 
that  soil  cover  does 
attenuate  overpressures 
appreciably;  mainly  due  to 
an  arching  action  in  the 
soil  as  the  membrane  yields. 
Had  these  models  been 
rigid  *hey  would  have  been 
subjected  to  pressures  dose 


Figure  18.  Square  Membrane  Test 
Yielding 


Figure  2C  .  Blast  Simulator 


Figure  19*  Square  Membrane  with 
Two  Yielding  Beams 


Figure  21 .  Model  Structure  and 


to  the  surface  overpressure. 


Deformed  Model  Shells 


As  con  be  seen  from  the 
test  evidence,  the  arching 
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ocficn  was  acting  well  before  the  model  was 
buried  one-holt  the  diameter  and  at  one 
diameter  only  one-tenth  of  the  overpressure 
is  felt  by  the  buried  membrane  roof  (Figure  22). 

These  results  are  not  useable  for  pro¬ 
totype  predictions  because  the  principles  of 
similitude  are  not  satisfied.  They  do  serve 
to  illustrate  the  attenuation  due  to  soil 
arching.  The  arching,  in  this  ccse,  has.as  a 
point  of  support,  the  edge  ring  of  the  model 
structure.  Had  this  edge  ring  not  been  °- 

present#ihe  soil  would  have  arched  from  a 
point  outside  the  structure;  but  the  arch 
would  have  been  longer  and  flatter  and  thus 
the  attenuation  of  overpressure  would  have 
been  less. 

In  tests  with  the  circular  models,  it 
was  noted  that  the  deflected  shape  was  nearly 
spherical  except  that  near  the  edges  where 
the  soil  arch  was  supported  by  the  structure 
the  curvature  was  greater.  The  rac’ii  of 
curvature  along  a  diameter  was  measured 
and  by  making  use  of  the  basic  formula 

1  i  _  P 

r  +  c  ~~~s 

I  z 


a  fairly  accurate  distribution  of  the  pressure 
across  the  membrane  was  determined 
(Figure  23).  Note:  a  thin  rubber  sheet 
was  placed  over  the  surface  of  the  sand  to 
prevent  the  blast  wave  from  permeating 
the  pores. 

Tests  have  shown  that  yielding  membranes 
have  the  ability  to  deform  dynamic  overpressures. 
It  has  been  noted  in  the  test  conducted  at  the 
University  of  Arizona  that  there  is  10  to  20 
percent  increase  in  deflection  under  a  dynamic 
load  as  compared  to  the  same  magnitude  of 
overpressure  applied  statically. 


Figure  23  Radius  of  Curvature  and  Pressure  Distribution  across 
a  Circular  Membrane  with  Soil  Cover 
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CHAPTER  IV 


SUMMARY  AND  CONCLUSIONS 


Yielding  membranes  as  elements  of  shelters  to  protect  civilian  population  from  overpressure  effects  of  nuclear 
weapons  are  definitely  feasible.  As  indicated  on  page  8  of  this  report,  under  certain  conditions  a  1/4  in.  thick  steel 
membrane  will  carry  more  load  than  a  19  in.  thick  concrete  slab  which  is  reinforced  with  steel  at  a  rate  of  4.08 
The  membrane  itself  uses  less  steel  than  the  reinforced  concre'e  slob.  Favorable  soil-structure  interaction  effects  are 
induced,  by  these  elements,  in  all  situations. 

Construction  difficulties,  brought  about  by  pre-dishing,  may  be  overcome  ns  more  information  is  obtained  on 
the  response  of  flat  membranes  to  blast  overpressures.  Problems  associated  with  corrosion,  and  continuity  of  welds,  may 
be  largely  overcome  by  using  (as  pseudo-membranes)  thin  concrete  slabs  which  are  reinforced  by  a  closely  spaced  wire 
mesh  which  runs  continuously  in  both  directions.  The  wire  mesh  then  acts  os  the  membrane  and  the  concrete  as  the 
corrosion  resisting  and  local  transfer  medium. 

Boundary  supports  and  full  scale  testing  are  items  that  need  more  attention  and  ore  suggested  as  parts  of  future 
investigations.  In  addition,  analytical  studies  on  strain  variations  and  time  response  functions  to  dynamic  loads  are 
needed . 


At  the  present  time,  however,  the  response  of  these  structures  has  been  sufficiently  bounded  to  warrant  their 
use  in  civil  defense  situations.  The  methods  used  in  this  feasibility  study  are  as  valid  end  reliable  as  any  that  are 
currently  being  used  on  the  design  and  analysis  of  conventional  types  of  structures.  It  is  recommended  that  the  Office 
of  Civil  Defense  seriously  consider  adding  these  types  of  solutions  to  the  ever  expanding  inventory  of  such  solutions* 
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APPENDIX  A 


YIELDING  MEMBRANE  FORCES 


Introduction 

This  Appendix  considers  the  configuration  taken  by  yielding  membranes  under  uniform  lateral  pressure  loading. 

It  deals  with  rectangular  membranes  which  are  rigidly  supported  on  four  sides  and  rectangular  membranes  which  have  one 
or  more  sides  supported  by  yielding  supports,  in  the  form  of  beams  which  deflect  as  the  membrane  deflects.  Some  of  the 
results  of  this  effort  were  presented  in  Chapter  III,  but  the  complete  details  ore  included  in  this  section. 

The  deflections  of  the  membranes  are  very  large,  producing  center  deflection-to-span  ratios  as  large  as  0.20 
which  results  in  average  biaxicl  strains  as  large  os  10  percent.  The  result  is  that  bending  forces  in  the  membrane  are  of 
little  importance  and  membrane  forces  predominate. 

If  D  is  defined  as  the  short-span  distance,  the  types  of  membranes  considered  are:  (1)  D  x  D  with  four  edges 
clamped,  (2)  1 .5  D  x  D  with  four  edges  damped,  (3)  2  D  x  D  with  four  edges  clamped,  (■*»,  D  x  D  with  a  plastic  edge 
beam  on  one  side  and  clamped  on  the  other  three  sides,  (5)  D  x  D  with  plastic  edge  beams  on  two  sides  and  clamped  on 
two  sides,  (6)  D  x  D  with  plastic  edge  beams  on  three  sides  and  one  side  clamped,  and  (7)  plastic  edge  beams  on  all  four 
sides.  Plan  views  of  the  shapes  are  shown  in  Figure  A-l . 

Only  membrane  forces  are  considered.  Further,  the  assumption  is  made  that  the  same  membrane  stress  level 
exists  at  all  points  on  the  resulting  surface.  The  problem  considered  in  this  investigation  involves  uniform  lateral  pressure 
only.  However,  the  differential  equations  developed  can  be  solved  for  other  types  of  lateral  loadings. 

Yielding  Membrane  Theory 

Formulation  of  the  Problem.  To  oroperly  predict  the  configuration  of  yie'ding  membranes  under  normal  pressure 
loading,  it  is  necessary  to  develop  adequate  theories.  Few  of  the  theories  presented  in  Historical  Review-  (Appendix  D) 
are  directly  applicable  to  the  problem  at  hand.  However,  miror  modifications  to  them  allow  appliability ,  The  work 
which  follows  presents  these  modifications. 


NOTATION 

c .  Width  of  plote  or  radius  of  circular  plofe 

b .  Length  of  plate 

c,d .  Dimensions  of  membrane 

D  .  Spon  of  membrane  or  flexuroi  rigidity  of  plate 

E  .  Modulus  of  Elasticity 

£  .  Strain 

F  .  Force  in  yielding  beam  or  stress  function 

h .  Thickness  of  plate,  x  direction  grid  size 

i  .  Grid  poirt  in  x  direction 

i  .  Grid  point  in  y  direction 

k .  y  direction  grid  size 

L .  Length  of  membrane 

N  .  Normal  force  per  unit  length 

P .  Pressure 


Poisson's  ratio 
Radius  of  curvature 
Membrane  strength 


r 

R 

S  . 


A-l 


<r- 

t  . 


Stress 

Constant  stress  (usually  the  yield  stress) 
Thickness  of  membrane 


I 

I 


u,vfw....  Components  of  displacements  in  x,y,z  directions,  respectively 
x, y , z  . . .  .  Rectangular  coordinates 

.  Vertical  deflection  in  center  of  membrane 


(- — D  — f“ - 1.5D - *|  \“  2D 


hH 

Figure  A-1 .  Plan  Views  of  the  Shapes  Considered 
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The  tangent  to  a  curve  at  a  point  M  is  the  limiting  position  of  the  secant  through  M  and  a  point  Mj  of  the 
curve  M)  approaches  M  as  a  limit.  The  equations  of  the  tangents  through  M  und  are 


dx  dy  dz 

Ti  3T  3T 


where  x,  y,  and  z  are  coordinates  of  a  point  on  the  curve,  their  values  depending  upon  the  parameter  h. 

*  H.  P.  Harrenstien,  "Configuration  of  Shell  Structures  for  Optimum  Stresses,"  Proceedings  of  the  Symposium  on 
Shell  Research,  Delft,  The  Netherlands,  196). 
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A  tangent  line  to  a  curve  upon  o  surface  is  called  a  "tangent  line  to  the  surface"  at  the  point  of  contact.  It  is 
evident  that  there  is  an  infinite  number  of  tangent  lines  to  a  surface  at  any  point.  However,  all  of  these  lines  lie  in  a 
plane  called  the  "tangent  plane"  which  is  tangent  to  the  surface  ct  the  point.  If  the  equation  of  the  curve,  C,  in 
curvilinear  coordinates  is  v  =  f(u),  then  the  above  equations  may  be  written 


x  -  x  =  h 
o 

d  X 
du 

+  '  ~r 

du 

■a r 

(6) 

y  -  y0  s  h 

6  y 

d  u 

>  f^y 

d  v 

du 

“3 r 

(7) 

Z  -  Z  -  h  I 

0 

1  ^  Z 
[  d  u 

*  f'  ±i 

d  v 

1  du 

I-*- 

(8) 

where  the  prime  indicates  differentiation.  In  order  to  obtain  the  locus  of  theif 

i  tangent  lines,  eliminate  f1  and  h  from 

these  equations 

.  This  gives 

X  -  X 

0 

y-y0 

z  -  z 
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d  x 

d  y 

d  z 
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6  u 
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d  z 

Tv- 

ST 

ST 

which  is  evidently  the  equation  of  a 

plane  through  M. 

The  equation  of  the  tangent  plane  may  be  written  as 

(x  -  X  )  X  - 
0 

(y  -  y, 

)  Y  *•  (z  -  ; 

0 

z  )  Z  =  0 

0 

(10) 
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The  positive  direction  of  the  normal-to-the-tangent  plane  is  defined  to  be  that  for  which  the  functions  X,  Y  and  Z  are 
the  direction  cosines. 


Consider  any  curve,  C,  on  a  surface,  S,  through  a  point,  M.  The  direction  of  its  tangent,  MT,  is  determined 
by  a  value,  dv/du.  Let  9  denote  the  angle  which  the  positive  direction  of  the  norma!  to  the  surface  makes  with  the 


positive  direction  of  the  principal 

normal  to  C  at  M, 

angles  being  measured  toward  the  positive  binormal  „  Thus 

Cos  9 

r  j  X 

d2x 
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Y  d2y 

Z 
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where  r  is  the  radius  of  curvature 
the  forms 

of  the  arc 

of  C.  In  terms  of  du/ds  and  dv/ds. 

the  derivatives 

in  rhe  parenthesis 
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Now  define: 
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Making  these  changes 


Cos  Q 
r 


2  0 
D  du  -  2D1  d>j  dv  4  D"  dv 

E  da  +  2  F  du  dv  t  G  dv2 


05) 


Consider  the  tangent  curve  in  which  the  surface  is  cut  by  the  plane  determined  by  MT  and  the  normal  to  the 
surface  at  M,  callod  the  'Vioimol  section"  tangent  to  MT,  and  let  rn  denote  its  radius.  Since  the  right-hand  member 
of  the  above  equation  is  the  some  for  C  and  the  normal  section  tangent  to  it. 


9  =  ±  (16) 

r  r 

n 

where  e  is  +1  or  -1,  depending  on  whether  9  is  less  than  or  greater  than  a  right  angle.  The  radii  r  and  rn  are  positive. 

Now,  let  us  introduce  a  new  function  R  which  is  equal  to  rn  when  0  12  5  */2,  and  equal  to  -r  when 
*/2  5  0  1  i  ,  ond  coll  it  the  "radius  of  normal  curvature"  of  the  surface  for  the  given  direction  MT.  As  defined: 


2  ? 

D  du  +  2  D1  du  dv  +  D"  dv 

- - - - - 

E  du  +  2  F  du  dv  +  G  dv 


or  on  substituting  t  -  dv/du 


I 

IT 


D  *  2  D1  t  +  D“  t2 
— 

E  *  2  Fm  Gt 


(17) 


(18) 


To  obtain  the  values  of  t  for  which  R  is  a  maximum  or  minimum,  differentiate  this  expression  with  respect  to  t 
and  set  the  result  equal  to  zero.  This  gives 


0  *  (O'  +  D"t)  (E  4  2  Ft  4  Gt2)  -  (F  ♦  Gt)  (D  ♦  2  D't  4  D"t2) 


(19) 


This  is  a  quadratic  in  t.  It  can  be  shown  that  at  every  ordinary  point  of  o  surface  there  is  a  direction  for  which  the 
radius  of  normal  curvature  is  a  maximum  and  a  direction  for  which  it  is  a  minimum,  and  they  are  at  right  angles  to  one 
another. 


Thus  the  two  values  of  R  become 


I  O’  4  D"t  ,  I  . 

ir  =  . t~+  gt  and  it 

and  the  following  relations  hold  between  the  principle 

E  4  Ft  -  R  (D  -  D't)  0 
F  *  Gt  -  R  (O'  4  D"  r)  -  0 


D  4  D't 
E  4  ft 

radii  and  the  corresponding  values  of  t 


(20) 


(21) 


Eliminating  t  from  these  equations 

0  -  (D  D"  -  D'2)  R2  -  (ED"  •  GD  -  2  FD')  R  4  (EG  -  F?)  (22) 

The  roots  of  this  equation  are  the  principal  «odii.  These  principal  radii  ore  denoted  by  and  R9  and  thus: 


A-6 


(23) 


V  7 


ED"  ‘  GD  -  2  FD* 

7 


T^ 


DD"  -  O' 


H 


Now,  for  ths  case  of  a  surface  defined  by  z  -  f(x,y),  fhe  following  values  are  obtained: 
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and  thus 
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also 


con' toot 


Using  this  exact  expression  for  the  normal  curvatures  tfie  equation  cf  equilibrium  becomes: 


(I  •  q  )  r  -  2  pqs 


•  (I 

TT 


2, 

P  >  » 


-  P 


(26) 


where  P  is  the  uniform  lateral  pressure  Jenote-d  by  previously. 

This  equation,  which  governs  the  shape  of  a  membrane  that  is  loaded  by  pressures  large  enough  to  coust 
deformations  in  excess  of  the  "large"  deflection  theory,  generally  may  not  be  solved  by  meons  usually  employed  for 
equctions  of  this  type.  The  equation  is  distinctively  non-linear,  involving  products  ond  powers  of  portiol  derivatives. 
The  method  presented  here  consists  of  defining  f!«  partial  derivatives  in  the  equotion  in  'em*  of  finite  differences 
ond  iterating  the  resulting  equations  by  Hie  Ncwfon-Raphson  feenr-ique  to  oefoin  u  description  of  th«-  surface . 
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An  oyer-relaxotion  procedure  was  used  to  speed  the  iterative  process.  This  iteration  procedure  was  too  involved 
to  attempt  by  desk  calculator  so  a  high-speed  digital  computer  was  used.  The  procedure  was  programmed  in  basic  Fortran 
computer  language  for  the  IBM  709  computer  at  the  University  of  Arizona.  The  computer  program  is  presented  in 
Appendix  B. 


Yielding  Edge  Beams.  In  :  'tuations  where  the  boundaries  ore  .ixed,  z  =  0  on  the  boundary  is  a  sufficient 
condition  to  admit  a  unique  solution.  However,  in  cases  where  these  boundaries  yield,  as  in  the  case  of  an  edge  beam 
which  yields  normal  to  the  initial  plane  of  the  surface,  this  condition  is  not  valid. 


In  this  case  the  differential  equation  governing  the  behavior  of  this  plastic  beam  across  the  membrane  is 
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PW  =  0 


(27) 


■where  F  is  the  plastic  strength  of  beam  plus  membrane  strength  and  W  is  the  beam  width.  The  second  term  is  a  load  term 
giving  the  contribution  of  the  pull  of  the  membrane  to  the  shape  of  the  plastic  beam.  The  2  in  the  numerator  arises 
because  the  membrane  is  considered  symmetrical  about  the  beam. 

This  equation  assumes  that  the  coordinates  extend  to  the  edge  of  the  beam.  Thus  F  is  actually  the  strength  of 
the  beam  and  the  membrane  combined,  which  must  to  taken  into  account  with  wide  beams  .  However,  in  the  computer 
study,  the  coordinates  were  taken  from  the  center  of  the  beam,  the  beam  width  was  assumed  zero,  and  F  became  the 
strength  of  the  beam  alone.  This  assumption  results  in  less  than  two  percent  error  in  most  cases. 

Finite  Difference  Equations.  The  equations  presented  for  the  membrane  and  the  yielding  edge  beams  were 
written  in  finite  difference  form  in  order  to  reduce  the  differential  equations  to  more  easily  handled  algebraic  equations. 
The  computer  programs  were  written  using  the  finite  difference  forms  of  equations.  In  the  programs  involving  the  yie!Jing 
edge  beams  the  two  sets  of  equations  were  solved  simultaneously. 


The  membrane  equation  written  in  terms  of  finite  d' Terences  about  a  point  involves  nine  points  as  shown  in 
Figure  A-2.  In  this  case  there  is  no  need  to  write  a  special  finite  lifrerence  equation  to  handle  points  along  the 
boundaries  because  the  boundary  points  are  either  zero  along  the  clamped  edges  or  a  point  on  a  yielding  beam.  The 
points  on  the  yielding  beams  are  obtained  by  solutions  of  the  beam  equation. 

The  partial  derivatives  written  in  finite  difference 

form  are 
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In  finite  difference  form  the  membrane  equilibrium  equation  becomes 
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For  a  square  mesh,  h  =  k; 
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I f  k  —  (d/c)h,  i.e.,  if  there  are  an  equal  number  of  divisions  i  each  direction,  then 
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Note  that  c  is  the  x-dimension  of  the  membrane  and  d  is  the  y-dimer.sion . 
In  finite  difference  form  the  yielding  beam  equation  is 
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F  Strength  in  beam  arid  membrane 

S  Strength  in  membrane 
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k  Grid  size  in  y  diiection 


♦  PW  «  0 


z.  Deflection  at  poiiit  i 


Computer  Results,  The  membrane  equation  in  finite  difference  form  (Equation  29)  was  soived  by  the  computer 
usiny  a  modified  Newfon-Raph:on  technique.  Pressure -to-membrane  strength  ratios  cf  from  zero  to  three  were  used  in 
the  computer  programs.  In  the  cases  of  membranes  supported  by  yielding  edge  beams,  Equations  29  and  32  were  solved 
simultaneously.  Values  of  the  ratio  of  beam  strength  to  the  product  of  the  membrane  strength  and  span  of  0.5,  1 .0  and 
2.0  were  used  in  the  programs. 

The  illustrations  in  this  section  are  the  result  of  the  computer  study.  Curves  are  computed  for  a  uniform  stress 
distribution  over  the  entire  surface.  However,  this  assumption  does  not  noticeably  affect  the  center  deflection  values 
as  will  be  shown  in  the  comparison  of  resuits  section.  In  fact,  the  only  creas  this  assumption  affects,  to  any  degree, 
is  the  comers. 


These  curves  of  PD/S  versus  z^D  are  general  curves  and  may  be  used  for  any  ductile  material  if  the  biaxial 
stress-strain  curve  for  the  material  is  available.  Only  in  the  case  of  a  rigid-plastic  material  can  the  S  value  be 
considered  a  constant. 


Figure  A-3  shows  the  standard 
grid  layout  used  in  all  of  the  following 
graphs,  tables  and  discussions.  It  was 
found  that  a  16  x  16  grid  gave  results 
within  the  desired  degree  of  accuracy. 

In  the  following  zj  .  means  zx  ^ .  The 
PD/S  versus  Zj/D  curves  (Figures  A-4 
through  A-9)  are  non-dimensionalized 
pressure  versus  center  deflection  curves 
for  membranes  with  clamped  edges.  The 
PD/S  versus  z/D  curves  (Figures  A-10 
through  A-16)  are  iKn-dimensiona!ized 
pressure  versus  deflection  curves  for  the 
points  of  maximum  deflection  of  the 
membrane  and  center  points  on  the  beams 
for  various  conditions  of  edge  constraint. 
The  numbers  inside  the  circles  are  values 
of  the  ratio  F/SD,  where 

F  =  Strength  of  beam  and 
membrane 

S  -  MenJbrane  strength 

D  =  Short  span  distance 

The  subscripts  on  the  circles  refer  to  the 
location  of  the  point  (see  Figure  A-3). 

For  example,  (T)o  9  refers  to  the  PD/S 
versus  z/D  curve  for  the  center  point  of 
the  membrane,  when  the  ratio  of  edge 
beam  strength  to  the  nroduct  of  the 
membrane  strength  and  the  span  is  one. 

The  graphs  and  tables,  shown 
relating  fa  membranes  supported  by 
various  yielding  edge  beam  arrangements, 
are  for  symmetrical  cases.  That  is,  it  is 
assumed  that  the  conditions  on  both 
sides  of  the  yielding  beams  are  the  same. 
These  tables  and  graphs  are  useful  in 
preliminary  design.  The  actual  values 
of  membrane  strength  and  beam  strengths 
involved  in  a  design  should  be  fed  into 
a  computer  program  to  determine  the 
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Figure  A-3.  Standard  Grid  Layout 
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Figure  A-4. 


PD/S  versus  z ^/D  Curve  for  a  Circular  Membrane 


TABLE  I 

Summary  of  Computer  Results;  Membrones  with  Clomped  Edges 
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Figure  A-9„  PD/S  versus  z  /D  Curve  for  on  oo  x  1  Membrane 
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Figure  A-l  1 .  PD/S  versus  z/D  Curves  fur  a  Square  Membrane  with  One  Yielding  Beam 


Figure  A-1A  PD  S  versus  z,  D  Curves  for  a  Square  Membrane  *ifh  Beams  on  Two  Opposite  Edges 
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TABLE  4  TABLE  5 

Summary  of  Computer  Results  -  Square  Membrane  with  Summary  of  Computet  Results  -  Square  Membrane  with 

Yielding  Beams  on  Two  Adjacent  Edges  Yielding  Beam*  on  Three  Edges 


Figure  A- 15.  PD/S  versus  z/D  Curves  for  a  Square  Membrane  with  Beams  on  Three  Edges  (Beams) 
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Figure  A-16.  PD  S  versus  z  D  Curves  a  <quure  Mcmtsrane  •■.ifr  Beams  or  Foi  r  £  x>e- 
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octual  deflections  of  the  surface. 
However,  olthough  the  edge  conditions 
and  strengths  and  loadings  of  one  panel 
affect  the  adjacent  panels,  it  not 
likely  that  they  will  affect  panels 
farther  away  to  any  significant  extent. 
Conservative  use  of  the  information  given 
here  will  allow  membrane  structures  to  be 
designed  within  the  desired  degree  of 
accuracy . 


Equivalent  Curves  for  Other 
Materials.  The  PD/S  versus 
curves  ore  completely  general  and  noy 
be  used  with  any  type  of  ductile 
material .  For  a  rigid-plastic  materiel 
the  S  value  is  a  constant  which  simplifies 
the  use  of  the  graph,.  If  these  curves 
are  used  a  number  of  times  with  a  certain 
material,  it  may  be  advantageous  to 
choose  some  constant  S  value  (such  as 
the  yield  strength)  and  draw  a  new 
"equivalent  rigid-plastic"  design  curve. 
This  eliminates  unnecessary  reference 
back  to  the  biaxial  stress-strain  curve, 
for  the  material,  to  obtain  an  S  value  at 
various  strain  love';. 

To  construct  this  "equivalent 
rigid-plastic"  curve,  a  stress  .alue,  0£, 
is  chosen  at  which  the  general  PD/  S 
versus  z^D  curve  is  to  coincide  with 
that  of  the  given  material  (pcint  A  of 
Figure  A- 17).  A  table  with  rhe 
following  headings  is  then  constructed 
(Table  7).  Columns  1,  2, and  3  ore 
obtained  by  '  Soos>  iq  a  number  of  points 
from  the  general  *  versus  zc/D  versus 
strain  curve.  A.  -..uji  three  or  four  points 
should  be  used  in  rne  portion  below 
PD/S  =  0.25.  From  the  biaxiai  »tresi- 
strain  curve,  stress  values  re  found 
which  coincide  with  the  strains  in 
column  3.  The  values  of  column  1  are 
multiplied  by  the  values  of  column  4 
and  divided  by  the  "constant"  stress  a^. 
A  plot  of  column  2  versus  column  5 
will  result  in  an  "equivalent  rigH- 
plasfic"  PD/S  versus  Zj/D  curve  for 
the  given  mcterial.  The  new  constant 
S  value  is  CFcf  • 


TABLE  6 

Summary  of  Computer  Results  -  Square  Membrane  with 
Yielding  Beams  on  Four  Edges 
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TABLE  7 

Column  Headings  for  " Equivolent-Rigid~Plas>ic’  Curves 
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Experimental  Analysis 

Introduction.  Tests  were  conducted  on  various  membrane  shapes.  These  tests  were  designed  to  obtain  membrane 
deflection  data  which  could  be  used  to  check  the  accuracy  of  the  computer  results  and  compare  with  the  predicted 
deflections.  Tests  were  run  on  1  x  i,  1.5  x  1  and  2  x  i  clamped  edge  membranes.  Tests  on  membranes  with  yielding 
beams  included  a  1  x  1  membrane  with  2  beams,  1 .5  x  1  membrane  with  2  and  3  beams  and  2  x  1  membranes  with  1,  3 
and  4  beams. 

Test  Equipment.  A  testing  device  was  constructed  to  test  the  membranes.  The  test  bed  was  a  25  x  15  x  3/4  in. 
high  strength  aluminum  plate.  Edge  beams  2-1/2  x  3/4  in.  thick  were  bolted  to  the  test  bed  by  1/2  in.  cap  screws 
spaced  1  in.  apart.  These  edge  beams  were  used  to  clamp  the  membrane  in  place.  Fluid  pressure  was  applied  to  the 
membrane  through  a  hole  in  the  test  bed.  One  edge  beam  could  be  moved  so  that  membranes  of  10  x  10  in.,  10  x  15  in., 
and  10  x  20  in,  could  be  tested  (see  Chapter  III). 

Two  tracks  ran  parallel  to  the  long  side  the  test  bed.  A  framework  was  fabricated  to  span  from  one  side  to  the 
other  over  the  membrane.  This  framework  held  a  sliding  bar  to  which  was  fastened  a  2  in.  dial  gage  which  measured 
within  0.001  in. 

The  dial  gage  was  positioned  crosswise  on  the  membrane  by  means  of  a  calibrated  screw  thread  which  moved 
the  sliding  bar.  The  entire  framework  was  positioned  lengthwise  on  the  membrane  by  means  of  another  calibrated  screw 
thread.  Great  care  was  taken  to  ensure  that  the  tracks  were  level,  the  sliding  bar  was  parallel  to  the  test  bed,  and 
the  gage  was  dead  centei  on  the  test  bed  when  the  screw  threads  were  positioned.  This  arrangement  permitted  test 
results  accurate  within  0.002  in.  in  the  vertical  direction,  and  within  0.005  in.  in  the  horteontal  direction. 

Water  pressure  was  applied  to  the  membrane  as  the  loading  medium.  This  reduced  the  danger  to  personnel 
which  could  have  been  present  with  air  pressure  only.  Air  under  pressure  stores  energy  which  would  have  been  suddenly 
released  when  the  shell  burst  and  could  have  caused  some  shell  fragments  to  be  thrown  about.  A  surplus  oxygen  tank 
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was  used  as  a  reservoir  and  air  pressure  to  water  pressure  conversion  unit.  A  sensitive  pressure  test  gage  was  attached  to 
the  tank.  Also  o  60  psi  pressure  celi,  in  conjunction  with  a  strain  box,  was  used  to  record  the  pressure.  The  pressure 
was  measured  to  within  0.015  psi. 

Test  Procedure.  Test  membranes  were  held  flat  by  a  styrofoam  block  as  the  edge  beams  were  attached.  The 
water  level  in  the  reservoir  tank  was  at  the  level  of  ihe  test  bed  at  the  start  of  each  test.  At  this  point  all  gages  and 
screw  threads  were  checked  to  be  certain  they  were  zeroed.  Deflection  measurements  were  taken  at  set  grid  points. 
(Each  membrane  was  divided  into  16  grids  each  direction  and  measurements  were  made  on  1/4  or  1/2  of  the  membrane 
depending  upon  the  symmetry  involved.)  As  a  result,  either  64  or  128  deflection  readings  were  taken  at  each  pressure 
level.  There  were  usually  five  or  six  pressure  levels  in  each  test.  Test  results  v^ere  reproduceable  and  so,  in  the  latter 
tests,  only  one  t  it  was  made  on  each  membrane  shape. 

Yielding  beams  were  simulated  by  No.  8  soft  steel  wire.  They  were  stretched  by  means  of  a  simple  jacking 
device  to  assure  that  they  were  well  into  the  plastic  range  before  the  test  was  started.  In  this  manner  the  load- 
displacement  curve  for  the  beam  was  essentially  horizontal,  and  the  load  in  the  beam  was  relatively  independent  of  the 
strain  and  could  be  more  closely  approximated.  All  the  initial  strain  during  testing  was  assumed  to  take  place  across 
the  width  ot  the  shell  because  the  edge  beams  clamped  down  on  the  wire  as  well  as  the  membrane.  This  strain  was 
obtained  from  the  Zq/D  versus  strain  curves.  At  points  where  the  wires  crossed,  a  small  device  bent  one  wire  over  the 
other  so  that  the  wires  themselves  were  in  the  same  plane.  The  added  strain,  caused  by  bending  the  one  wire,  was 
taken  into  account. 

The  yielding  beam  tests  were  all  of  a  type  in  which  sag  was  permitted  at  points  where  the  beams  crossed. 

The  apparatus  would  have  been  unnecessarily  cluttered  if  a  heavy  bar  had  been  used  to  hold  those  points  at  zero 
deflection. 

Experimental  Observations.  It  was  noted  that  as  the  membrane  deformed  the  material  seemed  to  stretch  quite 
unevenly  —  even  though,  at  the  same  time,  the  deflections  were  very  regular  and  predictable.  It  could  be  explained 
by  the  nature  of  the  thin  metal  membranes  --  copper,  aluminum,  and  steel  all  acted  the  same  way  —  in  that  the 
membranes  seemed  to  thin  out  along  irregular  paths  leaving  thicker  areas  surrounded  by  thinner  areas.  As  the  material 
was  stretched,  the  thicker  areas  gradually  decreased  in  size  and  finally,  just  before  failure,  most  of  the  earlier  thick 
areas  had  thinned.  It  seemed  that  this  thinning  of  the  materia!  occurred  at  first  where  the  strains  were  greater  (at  the 
center  and  about  half  way  toward  the  center  from  each  edge.  As  the  material  was  stretched  the  thinner  regions  were 
strained  into  the  strain-hardening  region  and  thus  reached  a  point  where  the  membrane  force  in  the  thinner  portions 
became  equal  to  the  force  in  the  thicker,  less  strained,  portions;  in  effect  equalizing  the  membrane  force  in  almost 
the  entire  membrane.  The  areas  in  the  corners  of  the  rectangular  shapes  were  the  last  to  be  strained  ir\  i  the  plastic 
region.  Also,  it  was  noted  that  the  circular  tests  (used  to  obtain  the  biaxial  stress-strain  curve)  the  straining  was 
regular  and  more  uniform. 

In  on  attempt  to  measure  the  strain  distribution  over  the  shell  surface,a  photo  grid  of  20  lines  to  the  inch  was 
fixed  on  the  surface  of  the  clamped  edge  membranes.  Near  the  edge,  the  strain  becomes  more  uniaxial  than  b'axial; 
therefore,  in  order  to  maintain  the  uniform  stress  level,  the  membrane  must  strain  more  near  the  edges  (see  Figure  A-18). 

As  is  shown  in  Figure  A-18,  the  stress  does  not  wry  much  even  though  the  strain  does  vary  considerably.  The 
points  in  the  square  symbol  are  the  stresses  taken  from  the  uniaxial  stress-strain  curve.  As  can  be  observed,  the  uniaxial 
stress  at  the  edge  is  the  same  value  as  the  biaxial  siress  at  the  center.  This  is  another  point  in  favor  of  the  uniform- 
stress-distribution  assumption.  The  significance  of  this  graph  is  also  considered  loter  with  material  on  comparison  of 
Computer  and  Test  Results. 

Because  of  the  irregular  strain  distribution  and  the  apparent  ability  of  the  material  to  equalize  the  membrane 
forces,  it  seems  that  the  assumption  of  a  uniform  stress  distribution  may  be  applicable.  This  assumption  obviously  results 
in  some  error  in  the  values  of  the  displacement  in  the  corners  of  a  rectangular  membrane.  The  amount  of  error  and  the 
region  it  involves  becomes  apparent  in  Computer  and  Test  Results  comparison. 

The  yielding  beams  strained  fairly  uniformly  across  the  span.  When  the  center  deflection  of  the  beam  become 
larger  than  one  inch,  straining  began  to  localize  near  the  edges.  Failure  did  occur  in  the  wire  during  one  of  the  tests 
and,  as  anticipated,  it  failed  near  the  edge  due  to  a  combination  of  tensile  and  bending  strains. 

The  stress-strain  curves  for  vhe  membrane  material  are  shown  in  Figures  A-19  and  A-20,  The  material  used 
was  annealed  steel  sheet  0.004  in.  thick.  The  load-strain  curve  for  the  yielding  beam  material  (No,  8  soft,  drawn. 
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Figure  A-19.  Uniaxial  Stress-Strain  Curve  for  Membrane  Test  Material 
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fencing  wire)  is  shown  in  Figure  A-21.  The  biaxial  stress-strain  curve  is  an  average  stress,  average  strain  curve 
computed  by  assuming  spherical  deformation  over  a  circular  membrane.  This  results  in  the  following  equations: 
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The  deformation  was  spherical  --  at  least  to  within  the  least  count  of  the  measuring  system  which  was  0,001  in. 
Although  the  strains  may  not  be  uniform,  they  were  assumed  to  be  for  the  computation  of  the  bioxial  stress-strain  curve 
because  this  was  to  be  the  assumption  used  in  the  rectangular  membranes. 

The  stress-strain  curves  and  the  membrane  tests  were  run  at  the  same  strain  rates  —  zero  strain  rate.  At  each 
load  level  the  deflections  were  allowed  to  stabilize  before  readings  were  taken.  In  this  way,  the  discrepancies,  which 
may  have  arisen  due  to  unequal  strain  rates,  were  minimized. 

The  deflection  readings  on  the  membranes  were  tkane  with  the  pressure  held  constant  so  the  deflections,  and 
thus  the  strains,  included  both  the  elastic  and  plastic  portions.  This  is  what  was  intended  since  a  normal  stress-strain 


Figure  A-20.  Biaxial  Stress-Strain  Curve  for  Membrane  Test  Material 
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curve  includes  both  portions.  However,  as  a  note  of  interest,  the  pressure  was  relieved  o  number  of  times  during  various 
tests  to  get  an  indication  as  to  the  effect  of  the  elastic  strains  on  the  deflections.  When  the  pressure  was  relieved  the 
decrease  in  the  amount  of  deflection  was  100  percent  up  to  z^D  -  0.023,  50  percent  at  =  0.030,  10  percent  at 
Zc/D  =  0.073,  5  percent  at  z&/D  ~  0. 103,  and  2  percent  at  z^D  ~  0.163.  These  can  also  be  obtained  by  using  the 
biaxicl  yield  strain  of  0.0C15  in/in  and  the  strain  versus  z^D  curves. 

Membrane  Boundary  Conditions.  The  membrane  forces  are  independent  of  bending  and  are  wholly  defined  by 
the  conditions  of  static  equilibrium.  However,  the  reactive  forces  and  deformation  obtained  by  the  use  of  the  membrane 
theory  at  the  shell's  boundary  usually  become  incompatible  with  the  actual  boundary  conditions. 

As  will  be  shown,  there  is  no  doubt  about  the  membrane  state  of  stress  existing  in  the  main  portion  of  the  shell. 
However,  very  near  the  boundary  the  clamped  edge  causes  a  narrow  band  in  which  the  curvature  is  of  the  opposite  sign 
from  that  of  the  main  portion  of  the  membrane.  In  this  region,  then,  there  are  large  bending  and  shearing  stresses. 
Fortunately,  this  region  is  of  the  order  of  thickness  of  the  membrane,  which  in  the  tests  was  only  0.004  in,;  or  more 
accurately,  the  boundary  layer  would  be^Rt",  which  in  the  tests  would  only  be  ^  ,03(,004)  in.^  ,or  0.01  in. 

Because,  in  this  edge  region,  bending  stresses  and  shear  stresses  exist  as  well  as  the  membrane  stresses;  it  is  very  likely 
that  this  is  where  failure  will  occur  first. 

If  the  strain  present  in  the  membrane  is  broken  into  two  parts,  due  to  tension  and  <2  due  to  bending,  then 
the  relative  amounts  of  each  can  be  determined  from  the  deflected  shape.  If  a  circular  deflected  shape  is  assumed  for 
this  purpose  it  will  be  a  simple  matter  to  determine  «]  and  <5*  Choosing  zc  -  center  deflection,  D  =  span, 

R  =  radius  of  curvature,  t  -  thickness  of  the  membrane,  then  it  can  be  shown  that 
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Figure  A-21.  Tensile  Load-Strain  Curve  for  Yielding  Beam  Test  Material 
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TABLE  8 

Strains  in  Center  Portion  of  Membrane 


For  the  tests  D  --  10  in.  and  t  =  0,004  in. 
Table  8  is  presented  for  the  center  portion 
of  the  membrane  shell. 


This  shows  the  insignificance  of  the 
bending  strain  <>ver  the  major  portion  of  the 
shell  area.  Also  if  Rw  R2  £  0  bending  stresses 
are  not  restricted  to  tne  edge  zone  of  the  shell . 
In  this  casef  however,  R  j  and  R2  are  both  of  the 
same  sign;  thus  Rj  Ro  >  0  and  the  bending 
stresses  are  restricted  to  the  edge  zone.  Just 
at  the  edge  of  the  membrane  the  curvatures 
measured  in  the  deformed  shells  were  as  small 
as  0.03  in.  This  occurred  at  the  midpoint  of 
the  longer  edge.  From  Figure  A-22 
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This  resuits  in  an  *2  at  the  center  of  the  clamped  edges  on  the  order  of  0.0625  in/in.  just  before  failure.  The  shells 
usually  failed  along  the  edge  when  the  z^D  ratio  was  somewhot  ever  0.20.  At  failure,  then,  the  total  strain  at  the  edge 
was  ‘  ]  +  *2  -  0.1635. 


This  rough  calculation  of  the  strain  existing  at  the  clamped  edge  at  failure  indicates  only  that  the  strain  would  be 
increased  by  the  flexure  occurring  at  the  boundary.  However,  using  a  Poisson's  ratio  of  1/2,  which  would  be  valid  in  the 
plastic  range,  it  can  be  shown  that  the  maximum  biaxial  strain  should  be  only  two-thirds  of  the  maximum  uniaxial  strain. 
This  explains  why  the  circular  membrane  tests  failed  at  strains  of  about  10  percent;  wheroas  the  uniaxial  test  specimens 
failed  at  strains  of  16  to  17  percent.  The  circular  membranes  did  not  fail  at  the  edges  but  rather  at  some  point  of  weakness 
in  the  material.  This  indicates  that,  although  the  rectangular  membranes  failed  at  the  edges,  they  were  near  failure  due 
to  biaxial  strains  in  the  center  portion  of  the  membrane. 


To  prevent  premature  failure  at  the  edges  of  a  membrane,  a  rounded  edge  of  some  type  must  be  provided,  h  the 
tests,  the  edge  beams  were  purposely  rounded  to  retard  failure  at  the  edges  and,  yet,  the  rounding  was  not  great  enough 
to  effect  the  deflection  readings  or  the  dimension  of  the  span. 


Analytical  and  Experimental  Comparisons 

Updating  Previous  Work.  From  published  results  of  ti  tj  conducted  on  circular  membranes  and  analytical  studies 
on  circular  membranes  an  accurate  check  of  the  PD/S  versus  z^D  curve  for  a  circular  membrane  may  be  obtained.  This  is 
shown  in  Figure  A-23  and  Tables  9  and  10,  Both  of  these  studies  attempted  to  acquire  stress  and  strain  distributions  over 
the  membrane  surface.  In  order  to  correlate  with  the  present  results  the  average  stress  in  the  membrane  was  used.  Note 
that  the  results  correspond  very  well  in  the  region  below  a  PD/S  volue  of  1,5;  obove  this  point  the  membranes  are 
approaching  the  "instability  strain"  or  the  strain  at  which  a  secondary  bulge  forms  at  the  crown.  The  greater  error  is 
found  in  this  upper  region  because  at  the  time  this  secondary  bulge  forms  a  noticeable  leveling  off  occurs  in  the  pressure 
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versus  deflecHon  curve.  The  PD/S  versus  z^D  curves  are  not  meant  to  account  for  this  loss  in  membrane  strength  ~ 
similar  to  the  decrease  in  stress  just  before  failure  on  a  stress-strain  curve  for  ductile  materials.  Failure  follows  soon  after 
this  "instability  strain"  is  reached. 

Test  results  for  the  deflection  of  square  clamped  plates  under  uniform  lateral  load  included  one  test  in  which  the 
center  deflection  to  plate  thickness  ratio  was  12:2.  The  tests  were  made  with  aluminum  specimens  with  the  following 
specifications:  D  =  7.5  in.,  t  =  0.0158  in.,  £  =  10,300  ksi,  CT  =  37.5  ksi  at  0.002  in/in  offset.  The  stress-strain  curve 
for  the  material  was  also  given.  Using  this  data,  the  5,  8,  10  and  12  w^h  ratios  were  converted  to  z^/'D  values  and  the 
Pa  /Et  ratios  were  converted  to  PD/S  values.  The  results  are  shown  in  Figure  A-24  and  Table  1 1 .  This  gives  some 
indication  of  the  point  at  which  the  bending  forces  can  be  taken  as  having  a  negligible  effect  on  the  deflections.  If  the 
bending  forces  were  affecting  the  load-carrying  capacity  of  the  plate,  the  points  plotted  should  be  above  the  curve.  The 
reason  the  points  are  below  the  curve  may  be  in  the  difference  in  strain  rates  between  the  stress-strain  curves  end  the  test. 
If  all  the  values  were  increased  so  that  the  largest  value  coincided  with  the  curve,  then  it  can  be  seen  that  somewhere 
between  the  Wg/h  values  of  five  and  eight  the  bending  forces  begin  to  exert  a  negligible  effect.  In  foct,  at  a  w^h  of 
5.1  the  bending  forces  carry  only  about  12  percenl  of  the  load. 

Results  of  failure  tests  on  1 .55  x  1  rectangular  membranes  with  clamped  edges  subjected  to  uniform  lateral  pressure 
are  shown  in  Figure  A-25  and  Table  12.  There  is  substantial  scatter  in  the  data.  The  ultimate  stresses  of  the  materials 
tested  were  obtained  by  uniaxial  tests  on  coupons  of  the  materials.  The  scatter  of  data  may  be  due  partly  to  the  difference 
in  the  rate  of  strain  between  the  coupon  tests  and  the  membrane  tests.  It  is  not  known  if  any  attempt  was  made  to  run  the 
tests  at  similar  rates  of  strain, 

Greenspon  considered  the  problem  of  large  deflections  of  a  plate  under  uniform  pressure  loading.  He  used  the 
Poisson  Equation  approximation 


TABLE  9 

Result*  of  Circulor  Diaphragm  Tests  (Gleyzal,  1948) 


Pa 

h 

(ksi) 

Ave.  Radiol 
Strain  (%) 

Ave.  Radial 
Stress  (ksi) 

W 

c 

a 

PD 

T 

Z 

C 

~u 

2.09 

0.035 

24.0 

0.024 

0.174 

0.012 

7.83 

0.190 

36.5 

0.054 

0.430 

0.027 

21.37 

0.875 

47.0 

0.113 

0.908 

0.057 

35.28 

1.580 

58.5 

0.156 

1.208 

0.078 

55.92 

3.200 

70.0 

0.221 

1  .595 

0.110 

71.68 

4.700 

74.5 

0.267 

1.925 

0.133 

84.76 

6.400 

78.5 

0.314 

2. 160 

0.157 

TABLE  10 

Results  of  Circular  Membrane  Tests  ('.Veil  and  Newmark ,  1955) 
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TABLE  11 

Results  of  Square  Membrane  Tests  (Romberg,  ef.  ol.  ,  19*12) 
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Figure  A-24.  Comparison  with  Ramberg,  et  ol  (1942),  Square  Membrane 


Figure  A-25.  Comparison  with  Greenspon  (1956);  Failure  Tests  on  1.55  to  l  Men^rane 
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Strain  (Percent) 


and  obtained  a  general  solution  to  the 
problem 

0.164  Po2 


TABLE  12 

Results  of  Failure  Tests  on  Rectangular  Membranes  (Greenspon,  1956) 


This  equation  may  be  rewritten,  in  term* 
used  here,  as  follows; 

Zc  0.164  PD 


The  comparison  of  this  Poisson  Equation  approximation  and  the  more  exact  equation  presented  here,  is  shown  in  Figure  A-26. 
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Comparison  of  Test  Results  and  Computer  Results,  All  of  the  comparisons  of  the  PD./S  versus  z/D  curves  made  in 
this  section  are  based  on  the  "equivalent  rigid-plastic"  curves-  Thus,  to  make  a  comparison  or  determine  the  magnitude 
of  the  errors  the  "equivalent  rigid-plastic1'  curve  (a  dotted  line)  and  the  test  points  (a  dot  in  a  circle)  should  be  used.  Use 
of  the  "equivalent  rigid-plastic"  curve  is  one  way  to  make  the  non-linear  behavior  of  the  material  more  apparent  (s»*e 
Figure  A-27  and  Table  13). 

It  will  be  noted  that,  for  all  the  shapes  tested,  the  points  in  a  Land  about  one-half  the  span  cf  the  membrane  in 
either  direction  across  the  center  of  the  membrane  show  negligible  error  between  the  test  results  and  computer  results.  In 
fact,  the  only  region  in  which  errors  appear  is  in  the  comers  in  an  area  with  a  rodius  of  about  one-fourth  the  span  with  the 
center  at  the  comer.  The  error  in  the  corners,  as  shown  in  Figure  A-28,  is  very  large  ot  low  pressure  levels,  and  about 
20  percent  at  high  pressure  levels,  This  is  to  be  expected  because  the  computer  program  assumes  a  uniform  stress  distribution 
which  does  not  exist  in  the  comers  of  u  rectangular  membrane.  However,  in  Figure  A-29,  for  a  point  o  little  further  away 
from  the  corners,  the  error  is  much  'ess  (less  than  10  percent  for  most  of  the  curve). 

In  Figure  A-30,  for  a  point  away  from  the  corners  but  along  the  edge,  the  errors  are  less  than  5  percent.  This 
same  trend  is  evident  in  oil  of  the  shapes  considered  and  also  in  the  membrane  supported  by  yielding  beams  (see  Figures  A-31 
through  A-40). 

In  the  illustrations  titled  ‘Comparison  of  Results  -  z/z^^..,"  beginning  with  Figure  A-41,  the  amount  of  error 
and  the  location  is  readily  apparent.  From  these  illustrations  a  few  characteristic  cross-sections  have  oeen  plotted. 

The  results  for  the  square  membrane  with  one  yielding  beam  show  excellent  correlation  between  computer  and 
test  results.  The  results  for  the  square  membrane  with  two  yielding  beams  show  an  error  of  about  eight  percent  in  the  main 
portion  of  the  membrane.  To  construct  these  charts,  first  the  stress-strain  curves  were  made  for  the  materiel.*,  used  in  the 
membrane  and  in  the  beams;  then  the  tests  were  conducted.  The  yielding  beams  were  stretched  to  about  three  percent 
strain  before  pressure  was  applied.  The  average  strain  across  the  point  of  maximum  deflection  was  computed.  Going  to  the 
biaxial  stress-strain  curve,  the  stress  in  the  material  was  found.  Thus,  the  values  of  pressure  and  membrane  strength  which 
were  fed  into  the  computer  program  were  determined  from  the  test  results.  In  this  way,  if  the  computer  program  gave 
results  very  nearly  the  same  as  the  test  results  the  contusion  could  be  made  that  the  assumptions  used  in  the  theory  were 
justified  (see  Figures  A-41  through  A- 53). 

It  can  be  observed  •  /,m  the  results  on  the  square  membrane  with  clamped  edges,  and  the  square  membrane  with 
one  yielding  beam,  that  the  results  of  the  tests  correspond  very  closely  with  those  of  the  computer  study.  The  strains  were 
easily  measured  in  both  of  these  cases,  so  accurate  values  of  the  membrane  strength  could  be  obtained.  The  errors  in  the 
square  membrane  with  tv  o  yielding  beams  were  due  mainly  to  the  errors  in  the  membrane  strengths  which  were  fed  into 
the  computer.  It  was  difficult  to  obtain  a  nigh  degree  of  accuracy  in  computing  the  strains  in  the  membrane,  and  thus  the 
membrane  strength.  The  program  is  very  sensitive  to  this  value  of  membrane  strength. 

The  greatest  errors  between  the  test  results  and  the  computer  results  occur  in  the  corners.  Here  the  membrane  is 
not  strained  as  much.  Because  the  strain  is  at  low  level,  the  stress  in  the  corner  is  also  at  a  lower  level  than  in  the 
main  portion  of  the  membrane,  and  in  fact  the  stresses  in  the  comers  approach  zero.  Thus,  because  the  strength  in  rhe 
comers  of  the  test  membrane  has  not  reached  the  same  strength  as  in  the  rest  of  the  membrone,  the  reflections  in  the 
corners  of  the  test  tnembrane  are  greater  than  die  computer  values  which  are  bused  on  a  ur  ioon  stress  distribution. 

As  can  be  suer  from  the  comparison  of  results  tobies,  the  two  sets  of  results  always  ogree  well  in  the  main  portion 
of  rhe  shell;  yet  have  errors  as  high  as  30  percent  in  the  comers  in  some  cases.  Note  olso  that  the  error  in  the  comers 
decreases  os  the  total  defection  increases.  This  is  because  the  stress  in  the  main  portion  of  the  membrane  has  been 
increasing  rapidly  in  the  clastic  ronge  of  the  stress-strein  curve  (see  Tables  14  and  15). 

Thu  membranes  yield  to  a  very  nearly  circular  deformation  pattern  in  the  short  spon  direction,  even  near  the 
clamped  edges.  The  greatest  deviarion  from  a  circulor  shape  occurs  at  o  distance  of  about  1  8  the  <pon  from  the  edges. 

A  circular  deformation  pattens  is  evident  from  the  strain  versus  center  deflection  curves  which  show  that  the  strain  varies 
very  little  f<M  the  different  rectongular  plans  and  the  circulor  plan.  The  circulor  membrane  and  the  strip  me  brone  both 
defiect  to  a  circulor  deformation  pattern  under  uniform  lateral  load. 

The  stress  and  strain  distribution  curve  (Figure  A-19hhcws  thoi  the  overog*  stress  is  about  44.8  Mi .  It  the  strain 
venws  z  0  curve  is  entered  with  i0  0  0.125,  the  value  of  strain  is  4.4  percent.  From  the  hio^iol  str-m-s train  curve 

(figure  A-10)  a  stress  value  of  48.5  ksi  is  obtained.  This  is  the  value  of  itress  used  to  obtoin  the  S  value  used  tor 


comparison  purposes  in  Figure  A-36.  Had  the  true  average  stress  value  of  44.8  ksi  been  used  the  correlation  would  have 
been  even  better* 


If  an  attempt  is  made  to  use  a  stress  function  which  has  a  maximum  value  ot  the  center  of  the  membrane  and 
approaches  zero  in  the  comers,  the  discrepancies  between  computer  results  and  test  results  will  decrease. 


Figure  A-27.  Comparison  of  Computer  ctnd  Tost  Results  -  Square  Membrane,  z, 
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*5,5 


A -39 


0 


0.04 


A  As. 


0  16 


0  12 

#/*> 

flgvn  A-4Q.  Co*npori«c*i  of  ComjHrtw  u-vd  T**t  Rwulh  -  2  *  1  Mtrrbto:* t 


A -40 


Sf-3  S-6* 
v  »  »  > 

bo  oo 

a .ps  r. « 

S  a*  oo  eft 

bo  '  b 

8*  8-S 


3  £  81  J-3  38 

bo  oo  bo  bb 

8-3  8-8  H  H 

bb  bb  oo  oo 

g. a  £.3  a  .a  s-a 

o  o  bb  bb  bb 

8-8  H  8-8  3-3 

bb  bb  bb  bb 


5  .8 


N 

oL 

■o  o 

S-3 

s-s 

U  »/> 

5-3 

(ft  is< 

s.  * 

w>  vS 

§*§ 

rs  go 

3*5 

8S 

8-8 

-or 

*-r 

b  o 

b  b 

b  o 

o  b 

o  o 

o  o 

o  o 

b  o 

3-5 

15 

m  to 

n  w) 

a*? 

8-8 

tx  O' 

a*  ft 

ss 

04  04 

5  .5 

PI 

i 

* 

-J 

n 

o  o 

b  b 

b  b 

b  o 

b  o 

b  b 

b  o 

b  o 

l 

Ol 

o 

04  04 

3.R 

CM  04 

fc.3 

04  04 

3.  a 

04  04 

$.8 

04  04 

3*2 

—  n 

W) 

s*8 

04 

K 

D 

; 

#» 

1 

o  b 

o  b 

b  o 

o  d 

b  o 

b  b 

o  o 

o  b 

» 

51 

Test 

u 

Jr 

X  O' 

GO 

rs. 

-0 

w> 

r> 

04  *- 

UUIOJ  Of  JO 


a. a 

o*  > 


a. ^ 

o  o> 


i  I  If3 

t~  U  M*  I 


1-2 
o  b 

a. 3 

O'  O' 

b  o 

2-S  a. 5f 

Os  Qk  CO  <30 

dodo 

S'? 

rs  >s 

o  o 

s-3 

b  o 

8-8 
b  b 

8-35 

04  04 

o  o 

0.892 

0.387 

■o  >o 

N  .  X 

00  CO 

o  o 

3-3 

o  b 

8.R 

Os  K 

o  o 

oo  rs. 

•  o 

K  'O 

o  o 

3*a 

o  o 

3*5 

b  b 

D  - 
fs. 

04  04 

o  o 

0.790 

0.788 

8*8 
b  b 

a. a 

o  o 

oo  rs 
—  •  O' 
x  « 

o  b 

rx  rs 

s/s 

b  o 

GO  04 
*5*04 
in  m 

o  o 

eg  O' 

3*s 

o  o 

04  04 

o  o 

5*3 

o  b 

8* 
o  b 

S-3 

o  o 

8-8 
b  b 

CD  Ox 
^  •  04 

m  «o 
d  o 

3*5 

o  b 

8-8 
b  b 

_  hs 

—  e  00 

cm  — 
d  o 

497 

483 

8- -3 

S-3  3-5 

eg  o> 

3-s 

2-8 

rs  — 
fx  ,vy 

CM  CN 

B-2 

i 

5 

o  o 

o 

o 

• 

o 

o 

b  b 

b  o 

•  • 

o  o 

0 

0 

0 

0 

*  * 

3 

fc.3 

8 

8. 

* 

rs  — 
rx.  * 

—  rx 

—  .00 

a. 5 

04  •  O 

1  " 

$  1 

o 

(N  (N 

o  o 

<N 

o 

04 

d 

04 

o 

CM 

o 

04  04 

o  d 

04  04 

o  o 

04  — 

o  o 

b  b 

o  o 

>o  to 

««!°d  P!J0 


A-41 


Figure  A-41.  Con*>arison  of  Results  -  z/i^^  for  Stluare  Membrane  Figure  A-42.  Comparison  of  Results  -  for  Square  Merriw 

with  Clamped  Edges,  z^y^  =  0.0594  with  Clamped  Edges,  z^yD  =  0.1265 
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Figure  A-43.  Cross-Sections  -  Square  Membrane  with  Clamped  Edges 
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FI9ure  A-45,  Comparison  of  Resulfs  -  z/z^  for  Square  Membrane  Figure  A-46.  Comparison  of  Results  -  z/z^  for  Square  Membran 

with  One  Yielding  Beam,  z^y'D  =  0.0575  with  One  Yielding  fleam™  /D  =  0.0937 
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Figure  A-49.  Comparison  of  Results  -  z/z  for  Square  Membrane  Figure  A-50„  Comparison  of  Results  -  Square  Membrane  with 
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with  One  Yielding  Beam,  z  /D  =  0 .1243  Yielding  Beams  on  Two  Adjacent  Edges,  z  /D  =  0.065 
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Figure  A-51.  Comparison  of  Results  -  Square  Membrane  with  Yielding  Figure  A-52.  Comparison  of  Results  -  Square  Membrane  with  Yielding 

Beams  on  Two  Adjacent  Edges,  =  0.1093  Beams  on  Two  Adjacent  Edges,  z  /D  =  0.1522 


APPENDIX  B 


YIELDING  MEMBRANE  PROGRAM 


Abstract 


This  Appendix  presents  a  computer  program  for  iterative  solution  of  the  large  deflection  membrane  equation 
using  finite  differences.  An  investigation  of  relaxation  factors  at  various  P^/S  ratios  and  grid  sizes  is  included.  Graphs 
ere  given  showing  the  relationship  between  strain  in  fie  membrane,  P3/S  ratios  and  center  deflections. 

Modifications  of  the  basic  program  are  given  which  allow  variation  of  pressure  across  the  entire  membrane,  use 
of  Poisson's  Equation,  inclusion  of  plastic  beams  across  the  membrane  along  the  center  lines,  and  creating  free  edges 
with  or  without  plastic  beams. 

Introduction 

The  equati  >n  defining  the  shape  a  membrane  takes  when  loaded  with  pressures  large  enough  to  cause  deformations, 
in  excess  of  those  allowed  by  small  deflection  theory,  defies  solution  by  means  usually  employed  for  equations  of  this 
type.  The  equation  is  distinctively  nonlinear,  involving  products  and  powers  of  partial  derivatives.  The  method  of 
solution  presented  here  consists  of  defining  the  partial  derivatives  in  the  equation  in  terms  of  finite  differences,  and 
iterating  the  resulting  equation  to  obtain  a  description  of  the  surface.  An  over-relaxation  procedure  is  used  to  spaed 
the  iterative  process.  The  iteration  procedure  involves  a  long  and  tedious  process  and  may  be  accomplished  only 
through  use  of  high  speed  digital  computers.  Therefore,  the  procedure  presented  here  is  Drogrammed  ir.  basic  Fortran 
computer  language. 


Mathematical  Formulation  of  the  Problem 

Consider  a  typical  element  of  surface  (Figure  B-l , 
£  j  and  are  general  curvilinear  coordinates,  which, 
when  combined  with  the  position  vector,  r,  describe 
the  surface  completely,  n  is  the  normal  to  the  surface 
at  the  point.  Then  the  general  equation  for  the 
principal  curvatures  kn  of  the  surface  as  developed 
by  Wang*  is: 

H2k2  -  (EN  -  2FM  +  GL)k  +  (LN  -  M2)  =  0  (1 

n  n 


El  *  constant 


f2  =cons 


Figure  B-l.  Typical  Element  of  Surface 


Solving  this  equation  for  k^  and  k^  by  the  quadratic  formula: 


3K  -  2FM  +  GL  +W(SN-  +  GL)  -  (IJJ  -  '/  ) 
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Tied  Elasticity,  McGraw-Hill,  New  York,  1953. 
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If  we  choose  for  our  generalized  coordinates  ^  =  x  and  ^  =  Yt  equations  (3)  become: 
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The  five  equations  developed  by  Wang  defining  the  equilibrium  of  the 
element  in  Figure  B-l  reduce  to  a  single  equation  if  ail  of  the  stress 
resultants  except  N;  and  N2  are  assumed  zero.  For  this  stress 
condition  to  exist  only  normal  loads  may  be  considered  as  shown 
in  Figure  B-2.  Let  and  R2  be  the  principal  racii  of  i-urvature, 
then  the  equation  of  equilibrium  of  the  element  i 


‘•'1  :h 

=0 


If  we  let  Nj  r-  N2  =  S  =  a  constant,  this  equaticn  becomes: 
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In  terms  of  principal  curvatures  this  becomes: 


nq 

^  k'j  ”  •  — *■ 

From  equation  (2)  we  see  that: 
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Substituting  equations  (4)  in  this  equation: 
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Figure  B-2.  Element  Under  Stress 


B-2 


Substituting  this  expression  into  equation  (5)  we  arrive  at  the  final  equations  describing  the  surface  of  a  membrane  under 
a  pressure  loading  p^. 
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The  finite  difference  expressions  for  the  partial  derivatives  in  equation  (6)  are: 
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Substituting  these  in  equation  (6)  and  reducing  we  find: 
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Where: 
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grid  sire  parallel  to  the  Y  «x  is 


If  equation  (8)  it  satisfied  by  the  correct  value  of  Zj,-  the  identity  will  hold.  However/  through  the  iterative 
process  the  correct  value  of  Zj  •  are  not  known  and  equation  (8)  is  not  equal  to  zero  but  rather  to  some  value  f  which 
approaches  zero  as  the  values  6f  Z*^  approach  their  correct  values. 

This  suggests  the  following  relaxation  scheme: 


(new)  «  (old) 
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<J>  *  ovcr-rolaxocton  factor 

»  -  2(4  X2  A2  ♦  4  X  2  +  B2) 

*  y 

Discussion  of  the  Computer  Program 

Figure  B-3  defined  the  various  parameten 
used  in  this  program  to  describe  the  physical  ** 

dimensions  of  the  n«mhrane.  Because  the  pressure 
is  assumed  constant  over  the  surface/  symmetry  can 
be  used  and  only  cne-quarter  of  the  membrane  i 

need  be  computed .  ' 

A  series  of  trials  was  made  to  determine 
the  optimum  relaxation  factors  that  would  give 
convergence  in  the  least  number  of  iterations. 

The  trials  were  made  at  PD/S  ratios  of  one 
and  two  (Figures  B-4  and  B-5).  It  should  be 
noted  that  the  relaxation  factor  is  not 
critical  for  large  grid  spacing  and  low 
PD/S  values  but  becomes  increasingly 
critical  as  either  PD/S  rises  or  grid 
spacing  decreases.  There  is  little  change 
in  the  relaxation  factor  with  respect  to 
the  convergence  limit  CONVL;  however, 
the  number  of  iterations  increase  at  the 
convergence  limit  decreases. 


i 


T  I 

Figure  B-3.  One-quarter  of  the  Membrane 


Figures  B-6  and  B-7  indicate  that  the  error  decreases  sharply  os  the  grid  size  decreases.  Less  than  one  percent 
error  may  be  achieved  with  an  8  x  8  grid.  These  curves  are  for  a  square  membrane  and  the  "exact"  values  shown  are 
from  a  32  x  32  grid  size  with  a  converg?'  <‘»  limit  of  CONVL  =  0.0001 .  This  was  considered  to  be  exact  enough  for  the 
purposes  of  an  error  analysis.  A  id  x  16  grid  size  with  CONVL  -  0.001  was  used  to  obtain  the  compc*er  PD/S  versus 
curvet.  It  was  found  that  minus  signs  appearing  in  the  output  column  giving  maximum  change  for  eoch  iteration 
indicate  the  relaxation  factor  is  too  high.  The  ideal  relaxation  foctor  was  found  to  produce  wily  one  or  two  minus 
signs  in  this  column. 


The  program*  hod  numerous  comment  statements  to  aid  in  identification  of  the  program  components  and  the 
variables  involved,  it  olio  hod  mony  additional  features  such  as;  initial  values  could  be  reod  in  os  data;  all  points 
could  begin  as  zero;  a  limit  on  the  number  of  iterations;  o  convergence  limit  (usually  CONVL  *  0.001  in.);  a  stop 
on  the  iterations  if  the  center  deflection  become  larger  than  one-half  the  span;  and  strain  computed  across  the 
point  of  maximum  deflection  ond  along  the  yielding  beams. 


Relaxation  Factor 


Grid  Size 

Figure  B-4.  Relaxation  Factor  versus  Grid  Size,  PD/S  =  1.0 
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Iterations  Iterations 


Capacity.  The  program  will 
handle  a  fini te-difference  grid  64  by  64 
on  one-quarter  of  the  membrcne .  The 
capacity  of  the  program  may  be 
altered  by  replacing  the  parameters 
of  the  DIMENSION  statement, 
appearing  at  the  beginning  of  the 
program,  with  new  parameters;  and 
replacing  the  limiting  values  of  the 
two  DO  statements,  immediately 
following  the  DIMENSION  statement, 
with  the  new  parameters.  The 
parameters  to  be  entered  must  be 
the  desired  maximum  size  of  grid 
plus  two.  This  requirement  arises 
because  the  program  includes  the 
boundary  points  and  one  line  of 
points  post  each  centerline.  Thus, 
the  altered  statements  would  read 
as  follows: 


DIMENSION  Z(XXX,XXX> 
C  ZERO  Z(I,  J)  BAND 

1  D02I  o  1.XXX 
D02J  -  l.XXX 

2  Z(I,  J)-0. 


where  the  new  array  size  would  appear 
in  piece  of  the  X's. 


Figure  B-6.  Percent  Error  versus  Grid  Size,  PD/S  =  1 .0 


Detailed  Procedure. 


Before  each  set  of  data  cards  is  read,  the  storage  spaces  for  Z(i,  J)  are  set  to  rero. 


Input:  Necessary  data  is  read  into  the  computer  on  two  data  cards. 

Card  1:  IDENT . .  Date  identification,  Kay  be  any  three 

dife,it  number.  It  will  be  printed  out 
with  the  results  as: 

IDENTIFICATION.  .  .XXX 

ICCOZ  .  This  parameter  tells  the  computer  whether 

or  not  initial  values  of  7.(1,  J)  are  to  be 
read  in  from  data  cards.  Punch  column  10, 
1  if  initial  Z(I,  J)  are  given,  2  if  not 
yiven. 
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Percent  Error  Percent  Error 


IX . The  number  of  finite  divisions  on  one- 

quarter  of  the  aenbrnne  parallel  to  DIMA. 
Max.  -  64 

IT . The  number  of  finite  divisions  on  one* 

quarter  of  the  membrane  parallel  to  DM. 
Max.  »  64 

ITERL . The  maximum  number  of  iterations  allowed. 

CGNVL . The  maximum  change  between  iterations  to  be 


allowed  at  any  point.  The  membrane  will 
converge  approximately  to  an  accuracy  one 
place  less  than  the  convergence  limit  set. 
Thus,  if  COHVL  -  .001  the  results  will  be 
accurate  to  .01. 


PHI  .  .  Over -relaxation  factor 

Card  2:  DIMA . Dimension  of  one-quarter  of  the  plate 

parallel  to  the  X  axis. 

DDE . Dimension  of  one-quarter  or  the  plate 

parallel  to  the  Y  axis. 

STRE3  .  The  membrane  stress.  The  stress  usually 


used  here  is  the  plastic  strength  of  the 
material  to  be  used  per  unit  length. 

Thus,  if  steel  with  a  yield  point  of 
36,000  psi  end  a  thickness  of  .125  inches 
Is  to  be  used ,  we  would  compute; 

STUBS  -  36,000  x  .125  -  4300  lbs/in. 

PRES3  ......  Pressure  acting  normal  to  the  membrane 

surface. 


It  is  possible  to  read  initial  values  of  deflections  into  the  computer  through  use  of  ICODZ.  (f  ICODZ  is 
punched  1  in  column  10  of  card  1,  the  computer  will  read  in  initial  values  for  deflections.  If  ICODZ  is  punched  2  in 
the  same  column,  the  ccmputer  will  bypass  reading  of  initial  deflection  values  and  'nitial  deflections  will  be  assumed 
zero.  If  initial  deflections  are  to  be  read  in,  the  follow:ng  procedure  must  be  followed.  The  machine  will  read  the 
first  word  of  the  first  card  as  Zjj,  the  second  as  Z^/  etc.,  until  the  first  row  has  been  read  in.  The  machine  finds 
the  nuinber  of  pieces  of  data  to  place  in  the  first  row  by  interpreting  the  parameters  given  it  on  Cards  1  and  2.  If  more 
than  one  card  is  required  to  fill  the  first  row,  the  machine  will  read  n  second  card  ond  continue  placing  numbers  in  the 
first  row  until  the  first  row  has  been  completed.  The  first  number  of  tine  second  row  will  then  be  read,  i  .e«,  ?-2y  as 
the  number  immediately  following  the  last  number  of  the  first  row,  A  new  card  must  not  be  started  for  each  new 
row  unless  the  first  value  for  that  row  begins  a  new  card  as  a  natural  port  of  the  sequence.  The  data  must  be  placed 
on  the  cards  in  the  10  column  fields,  with  decimal  points  included,  continuously  beginning  with  the  first  deflection 
of  the  first  row  ond  ending  with  the  center  point  of  the  membrane.  All  boundary  points  must  be  included  but  only 
one-quarter  of  the  membrane  m ay  be  read  in.  The  following  example  demonstrates  this  procedure.  Suppose  we  have 
the  plate,  given  in  Figure  8-8,  with  initial  deflections  as  shown  at  the  grid  points.  Then  suppose  we  wish  to  determine 
the  deflected  shape  of  the  plate  under  pressure  of  20  psi.  The  plastic  strength  of  the  plate  is  100  Ibs./in.  The 
information  would  be  entered  on  the  data  cards  as  follows: 

Card  1; 


ITEM 

CARD  COL. 

EXAMPLE 

IDENT 

8-10 

xx  xxxx  x!25 

ICCDZ 

20 

xx  xxxx  jr-xl 

IX 

29-30 

xx  xxxx  so. 04 

IY 

39-40 

xx  xxxx  xx  04 

ITERL 

41-50 

xx  xxxx  xx50 

C0NVL 

51-60 

soc  xxxx  .001 

PHI 

61-70 

sex  xxxx  1.50 
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Card  2: 


DIHA 

1-10 

XX 

xxsx 

xx5. 

DIKE 

li-20 

XX 

xxxx 

xx  5 . 

5TRES 

21-30 

XX 

xm 

100. 

TRES3 

31-40 

XX 

x20. 

CARD 

C<' RTS  CCt.. 

ronrr 

EXAMPLE 

3 

1-10 

2i.i 

00  0000  0000 

3 

11-20 

Z1.2 

00  0000  0000 

3 

21-30 

21.3 

00  0000  0000 

3 

31-40 

Z1.4 

00  0000  0000 

3 

41-50 

*1.5 

00  0000  0000 

3 

51-60 

22,1 

00  0000  000c 

3 

61-7C 

Z2,2 

00  0000  00.2 

3 

71-80 

Z2*3 

00  0000  00.3 

4 

1-10 

Z2.4 

00  0000  00.3 

4 

11-20 

Z2.5 

00  0000  00.4 

4 

21-30 

Z3.1 

00  0000  0000 

etc. 

Where  x's  are  shown  the  card  columns  nay  be  left 
blank.  Decimals  may  not  be  included  where  not  shown. 

Where  decimals  are  shown  the  data  may  be  placed  anywhere 
within  the  ten  column  field  designated,  provided  the  decimal 
is  supplied.  If  the  decimal  is  not  supplied,  it  will  be 
automatically  placed  after  the  last  column  in  the  field 
regardless  of  where  the  data  is  placed.  Preceding  zeros 
need  not  be  punched . 

After  reading  in  all  dota  cards,  the  machine  computes 
necessary  parameters  and  sets  up  checks  for  use  later  in  the 
program.  The  f* rst  value  computed  is  Z^2 •  Equations  (9) 
are  evaluated,  using  the  initial  values  given,  then 
substituted  into  equation  (8).  The  value  of  f,  thus  found, 
is  substituted  into  equation  (10)  and  a  new  value  of  Z2  2 
is  found.  This  value  is  checked  against  the  old  value 
of  Z2  2  t°  see  if  convergence  has  been  obtained.  The 
old  value  is  then  replaced  by  the  new  value  and  the  program 
moves  to  point  2,3  where  the  process  is  repeated. 

At  each  point,  the  new  value  is  subtracted  from  the  old  value  and  the  difference  is  compared  with  the  largest 
difference  yet  found.  If  the  new  difference  is  larger,  it  replaces  the  old  largest  difference.  At  the  completion  of  an 
interation  the  largest  difference  anywhere  in  the  plate,  together  with  its  coordinates,  is  printed  out.  If  this  difference 
is  smaller  than  the  convergence  limit  set,  the  program  goes  to  print;  if  not,  the  program  returns  to  2,2  and  the  process 
is  repeated. 

At  the  center-lines  the  program  makes  use  of  symmetry  to  obtain  values  of  points  outside  the  quarter-plate 
being  computed.  At  the  completion  of  each  iteration  the  number  of  iterations  completed  is  checked  against  the 
iteration  iimit  given  in  the  input  data.  If  the  limit  is  reached,  before  convergence  is  reached,  the  program  prints 
out  all  the  data  obtained  to  that  point  and  goes  on  to  read  tne  next  data  cards.  Each  point  is  computed  consecutively, 
row  by  row,  beginning  with  point  2,2  and  ending  with  the  center  point  of  the  plate. 


DIMA  =  5” 

-3- 

“  -V-  =  1.25" 

1-  =  4 


Figure  B-8.  Example 
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After  the  equations  have  converged  or  the  limit  set  on  iterations  has  been  reached,  the  program  computes  the 
average  strain  along  the  centerline  parallel  to  DIMB  and  along  the  centerline  parallel  to  DIMA.  To  find  this  strain  the 
computer  sums  the  chord  distances  between  grid  points  along  the  centerline  and  subtracts  the  plan  projection  of  the 
centerline  from  it.  This  gives  the  elongation.  To  find  the  average  percent  strain,  the  elongation  is  divided  by  the  plan 
projection  of  the  centerline  and  multiplied  by  TOO. 

It  was  found  thar  the  equations  diverge  after  the  average  slope  (computed  as  the  center  deflection  divided  by  the 
shortest  side)  exceeds  1:1.  The  computer  computes  the  slope  along  both  centerlines.  The  program  will  automatically 
stop  iterating,  when  this  slope  is  reached,  and  print  out  all  available  data  plus  the  following  statement: 

AVERAGE  SLOPE  EXCEEDS  1:1,  EQUATIONS  DIVERSE  BEYOND  THIS  POINT 

Output.  A  flow  diagram  and  a  print  of  the  computer  program  follows.  The  program  is  supplemented  with  comment 
statements  to  facilitate  coordination  with  the  flow  diagram.  The  first  four  cords  of  the  program  (the  fourth  card  is  blank) 

Flow  Diogram,  Funicular  Shell 
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Flow  Diagram,  Funicular  Shell 


Flow  Diagram,  Funicular  Shell 


Return  Switch=2 
Z. .  *  •” Z.  > 


F low  Diagram,  Funicular  Shell 


Compute  Ave*  Strain 
at  X  =  a  and  Y  =  b 


Print: 

Iterations  Completed 
Z  at  Center  Point  of  Shell 
Point  of  Max.  Change 
Max.  Change 


/Print: 

Data  Identification 
Relaxation  Factor 
Ave.  Strain  at  X  =  a 
Ave,  Strain  at  Y  =  b 

p3,  s,  pys 


V 

'x' 


a,  b 


! 


Print: 

Z.  .  for  Entire  Shell 


f Print: 

Iterations  Completed 
Z  at  Center  Point  of  Shell 
Point  of  Max.  Change 
Max, 


Change 


_ i _ 

Add  1  to  Iterations 

_ s 

A 

***  HANSEN  2R 

**  PROGRAM  IS  ITERATIVE 

♦  COMPILE  FORTRAN*EXECUTE  FORTRAN 

C  FUNICULAR  SHELL*  PRESSURE  CONSTANT 
C  RESTRAINING  BANDS  ACROSS  X=0  AND  X=2*DIMA 

DIMENSION  Z ( 66 * 66  ) 

C  ZERO  Z ( I  •  J )  DANC 

1  D02 1=1*66 
D02J=1 *66 


2  Z( I  * J) =0. 

C  MAX  IX  =  MAX  IY  =  64 

C  ICODZ. ...PUNCH  COL  10.  1  IF  Z(I.J)  GIVEN*  2  IF  NOT  GIVEN 

C  IX . . .NUMBER  OF  DIVISIONS  ALONG  DIMA 

C  IY . .NUMBER  OF  DIVISIONS  ALONG  DIMS 

C  ITERL... .MAXIMUM  ITERATIONS  ALLOWED 

C  CONVL... .CONVERGENCE  LIMIT. ..MAX  CHANGE  BETWEEN  ITERATIONS 

C  DIMA . HORIZONTAL  DIMENSION  OF  1/A  SHELL 

C  DIMB . VERTICAL  DIMENSION  OF  1/A  SHELL 

C  STRES... .STRESS  PER  LINEAR  DISTANCE 

C  PRES3... .PRESSURE  PER  UNIT  AREA 

C  PHI. .....RELAXATION  FACTOR. ..USE  1.7  FOR  AN  8  BY  8  GRID 

C  STRNB... .YIELD  FORCE  OF  BEAM*  POUNDS 

C  ITER... ..ITERATIONS  COMPLETED 

C  Z( I .J5...RIZE  AT  POINT  I , J 

C  T CON2 • • • • MAX  CHANGE  BETWEEN  THIS  ITER  AND  LAST 
C  IP.JP... .POINT  OF  MAX  CHANGE 

C  PS . PRESSURE  OVER  STRESS  RATIO 

C  STNBX... .AVERAGE  STRAIN  OF  RESTRAINING  BAND  (BEAM)  AT  X  =  0 

C  STNSX... .AVERAGE  STRAIN  OF  MIDDLE  LINE  AT  X=DIMA 

C  COMPUTED  AS  100  TIMES  (ARC  LENGTH  -  DIMB) /DIMB 

C  XLAM . MESH  SIZE  ALONG  DIMA 

C  YLAM ... 4 . MESH  SIZE  ALONG  DIMB 

C  INPUT 


READ  50C.IDENT, ICODZ. IX.IY, ITER'  .CONVL. PHI 

500  FORMAT (5110. 2 FlU. 0) 

READ  501 .DIMA. DIMB, STRES. PRESS .STRNB, WIDTH 

501  FORMAT ( 6F 1 0  •  C  ) 

PRINT  599 

599  FORMA T(54H1FUNICULAR  SHE LL , RES TRA I N I NG  BANDS  ACROSS  X  =  0.  X  =  2DIMA/ 
H8H  I  DEN  T  IF  I  CAT  I  ON.  ..  I  3  ) 

600  FORMAT ( 75H  ITER  CENTER  DEFLECTION  MAX  CHANGE  IN  SHELL 
1  MAX  CHANGE  IN  kEAM/) 

C  COMPUTE  ARRAY  SIZE 

M=  IX+1 
N  =  I  Y+  1 

C  CHECK  TO  SEE  IF  INITIAL  Zd.J)  ARC  GIVEN 

I F ( I  COD? -  1 >9 99,3,4 
3  READ  50?,(  (Zd.J). J=l,y.), 1  =  1  .N) 

502  FORMAT  ( 6 F 1 o . u ) 

A  X  1  =  I  X 

Y  I  =  I  Y 

xlam  =  d;ma/x  i 

YLAM=DIMB/Y I 
i  T  E  P  =  1 
10  1=2 
J  -  2 

T  CON?  =  0 
T  C  ONA  =  0 • 

ICSWHd 
I  !  W  C  H  -  1 
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GO  TO  100 

101  J=J+l 

IFIM-Jl U .12,100 
12  2(1 ,J+1 1=2(1 ,J-1  1 
GO  TO  100 
11  1=1+1 
J  =  2 

I F  ( N- 1  1999,15.100 
15  I SWCH=2 

2  (  1  +  1 , J 1 =2 ( 1-1. J) 

GO  TO  100 

102  J  =  J+ 1 

I F ( M- J ) 1 50.20,21 
202(1 .J+l 1=2(1 .J-l 1 

2  (  !  +  l  .  J ) =2 ( 1-1. J) 

GO  TO  100 

21  2(1+1,01=2(1-1. J1 

GO  TC  100 

C  CCMPUT  A  TUCN  PH A o L 

c  compute;  nfv:  zti.j) 

IOC  A  =  2 (  I ,J+1  1-2(1 .J-l ) 

3  =  2 ( 1+1 , J 1  -  2  <  I-l.J) 

C  =  2( l  +  l »J)-2.*2( 1 ,J)+2( I  —  1 , J 1 
D=Z( I . J+l 1-2. *2 ( I , J l+2( 1 , J-l ) 

£  =  Z( 1  +  1 , J+ 1  )  -2 ( 1-1, J+l) -2 (l  +  l, J-l 1+ 2(1-1 , J— 1 ) 

F=<4.*XLAM»«2  +  A*»2)I!C-.5*A*P*L+<4.»YLAM**2  +  LS»*2)*D 
G  =  PRc03/  (  STit£o«  2  .  *XL  AM*  YIAV)*(4.*XLAM**2*YLAM**2+XLAM**2*B**2 
1 +YEAY* *  2  * A**2 1  *  *  1 .5 
l-.  =  2  .  *  (  4.  *  X'w'Ai*  +  «;+A  *  *2+4.  *  Y'_AV*»2  +  d**2  ) 

T  2  (  I  ,  J  1  +  Pill  *  ( I-  +  G  )  /h 
PS  =  PPES3/5>T  (ti 

C  GET  UP  CONVERGENC E  CHECK 
T  COM  =  T-Z  (I  ,  J  1 

I F I  APEC { TOCN1 l-CONVL) 111,111,107 
107  ICSWH«2 

111  I F ( A3SF ( TC0N1 1  -AoSF  (  T  COM2  1  1  1  Of. .  1  06 . 1 08 
ICe  T COM2  =  T CON  1 
1P=  1 
JP  =  J 

106  2 (  1  ,  J  1  =  T 

C  FIND  RETURN  STATEMENT 

109  GC  TO  (101 , 102 1 . ISWCH 
C  COMPUTE  BEAM 

150  1=2 
J  =  1 

155  AB  =2 ( I , J  + 1 5 -2 ( I ,J) 

3B  =  2  t 1  +  1 . J 1 -2 ( I  -  1 , J 1 

C5=Z ( 1+1 » J 1 -2.*Z  t I ,J)+2( 1-1 »J) 

DB=  ( 4 .  «  Y'_AV *«  2  +  tiiJ * *2  I  **1,5 
Ft>  =  2  • *5  TRES* AS/5GR  TF(XLAM**2+AO**2 1 
G6  =  8 . *  S T  RN6 *  V_  AM«CU/Ob 
GBD=PRES3*W1DTH 
Hb  =  16.«STRMu*YLAM/L>b 
TB  =  Z  (  I  ,  J  1  +  (  FB+GtJ+GBB  1  /HB*PHI 
C  SET  UP  CONVERGENCE  CHECK 
T  CON3  =  T  B-2 (  I » J ) 

I r  ( AH SF  (  T COM3  l-CO.NVL  1151, 151, 152 
152  ICSWH=? 

151  IFtABSF  t  TC0N31-ABSF ( TC0N4) 1153.153,154 
154  T C0N4  =  T COM3 

I  PB=  I 
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JPB*  J 

153  Z  ( 1  » J 1 *TB 
1  =  1+1 

1 F ( N— I )200»156»155 
156  2 ( 1+1 t J 1 =2 ( 1-1 *J ) 

GO  TO  155 
CHECK  CONVERGENCE 
200  GO  TO  <203,2021 .ICSWH 
CHECK  ITERATIONS 
202  IFt ITERL-ITER1203. 204.204 

DATA  TO  BE  PRINTED  OUT  EACH  ITERATION 
204  PRINT610.ITER.Z(N.M) • IP.JP.TCON2  *  I PB » JPB . TC0N4 

610  FORMAT) I5.F18. 5.  1  OH  2<  12 >  1H.  12.  2H!=  F9.5. 

1  10H  2 (  12*  1H.  12.  2H )  =  F9.5) 

CHECK  AVERAGE  SLOPE 

IF ( 1. -2 (N .Ml /DIMS  129.29. 30 

29  PRINT6U 

611  FORMAT ( 64H  AVERAGE  SLOPE  EXCEEDS  1  TO  1 ♦  EQUATIONS  DIVERGE  PAST  TH 
IIS  POINT  ) 

GO  TO  2C3 
ITERATE  AGAIN 

30  I TER= I TER+ 1 
GC  TO  10 

COMPUTE  STRAIN  IN  BEAM  AT  X=0 
2C3  1=1 
J  =  1 

ARC=0. 

41  APC1 =SORTF ( YL A”* *2+ (Z<  I  +  1 » J 1 - 2  < 1.J1 1**21 
ARC=ARC+ARC1 

i  =  r+i 

IFIM-l  1999. 62. 41 

42  STNBX= (  ARC-DIM.B  1  /D  I  MB*  100. 

COMPUTE  STRAIN  IN  SHELL  AT  X=A 
1=1 

J=M 
A  F.  C  =  0  • 

?1  ARC1  =SORTF  (  Y'.AK**2  +  <  2  (  1  +  1 ,  J)  -Z  (  I  .  J  1  1  **2  1 
ARC=ARC+ARC 1 
1  =  1  +  1 

I  FIN- 1  1999.32.31 
32  S I  NBX=  ( A  P.  C  -  J  I  P.'S  1  / 1/ 1 M  c  *  1*0  0  » 

PR  1ST  61  0,  I  TER  .Z  <M..”  1  .  IP.  JP.  T  CON  2  .  IPB.JPB  .TC0N4 
op  I 'IT  59°  ,  I CEN  T 
PRINT  609.PKI 

609  FCP’-'AT  (  1  7H  RELAX.  F  AC  TOR  •  .  .  F10.  5  1 
P'lNT  612.V.S1MSX 

612  FORMAT  H6H  ..TRAIN  ALONG  J=IZ.4H  IS  F3.2) 

PRINT  691  .PRESS.;.  TN:‘.X 

631  FORMAT  I  1  th  ........  F  1C  .  5 .2PH  STRAIN  IN  BEAM  .  .  .F 1 0 . 5  1 

PRINT  602  .r.I-tES.STKNiJ 

602  FORMATI17H  SHELL  o  T  <ENL>  1  H.  .  c  1 0 . 5 . 2  OH  riL-'v  ST  P.L  N'jTH  .  .  .  .  F  1  fl .  5  1 

PRINT  603.PS./.10TH 

603  FORMAT  (17H  P/S... . .  •  •  F  I  P  .  5  .  2  OH  PlAV  vJ  •  OTM  . F10.5) 

pp  I  NT  606,  XLA.v.  DIVA 

604  FORMAT  (  >  Yj(  LAM'.'A  X  .........  r  I  0. 6  •  2  OH  A  ................  F 1  0 . 5  1 

PP  INT6  '•?  ,  Y'.AV.q  I'.';! 

605  FORMAT  (17H  LA.n>6  Y . ...Fl  0.5*2  PH  I. . .  F 1  C  .  3  ) 

n»!MT  606. IX. I Y 

606  FORMAT  (  1  7H  L'rt  1  u  o!ZE . 13.41.1  ,ivl3) 

PRINT  6C7 

6C7  FORMAT I26MCZ  uloIANCtS  FOR  1/4  SHELL/) 


.F1P.5) 


■  F 1  C  .  3  ) 


p  p  I  f  ‘  T  6 \‘ ,  l  (  1  ,  J  »  i  (  I  .  J  )  .  J  •  1  .  I  .  I  *  1  .  Pi  > 
mi  crp.-AI  I  f:H  2  (  !  s  .  1st «  I  2  ,.7H  1  =  TB •  5 .  Sh 
1  5H  /.I  I.’,  l  i  >  »  I  2  .  ?H  1  =  F  F.  .  5  .  f-t : 


/.(  I  3 . 1H,  1  ?  ,2H)  =F8,5  « 
Z I  ! 2 . 1 H , I  2 ,2H1 =F6.5 * 


are  instructions  for  the  computer  necessary  for  correct  handling  of  the  program  and  for  time  charges.  These  cards  will, 
of  course,  vary  with  the  individual  system.  A  print  of  the  sample  data  follows  the  program.  This  data  was  run  on  the 
computer  for  both  cases  wherein  initial  values  were  given  and  not  given.  It  is  interesting  to  note  that  the  same  number 
of  iterations  was  required  for  convergence  in  both  cases.  Closer  investigation  of  the  output  for  each  iteration  revealed 
that  initial  rise  of  the  membrane  is  very  rapid.  Indeed,  by  the  end  of  the  sixth  iteration,  the  center  deflection  for  both 
cases  was  almost  identical.  This  property  makes  the  method  even  more  desirable  because  initial  values  are  not  necessary 
in  order  to  reach  convergence  within  a  reasonable  computation  time. 

The  first  page  of  output  gives  storage  locations  of  the  parameters  used  in  the  program.  This  information  expedites 
modifications  in  the  program  because  all  statement  numbers  and  all  variable  names  are  given.  The  second  page  of  output 
gives  the  center  deflections  at  each  iteration  and  the  point  of  maximum  change  and  the  value  of  the  maximum  change 
for  that  iteration.  The  third  page  of  output  gives  ail  of  the  final  dota. 

Relaxation  Factors 


A  series  of  trials  was  made  to  determine  relaxation  factors  that  would  give  convergence  in  the  least  number  of 
iterations.  These  trials  were  all  made  on  a  10  x  10  inch  plate  with  a  yield  strength  of  100  lb$./in.  The  trials  were 
made  with  pressure  of  10  and  20  psi.  Figure  P-9  gives  a  comparison  of  the  two  pressures.  It  should  be  noted  that  the 
relaxation  factcr  is  not  critical  for  low  grid  sizes  and  low  pressures  but  becomes  increasingly  critical  as  either  the 
P-jD/S  ratio  rises  or  the  grid  size  increases.  For  a  16  x  16  grid  at  a  P^/S  ratio  of  2.0  the  graph  is  extremely  shape  and 
slight  error  in  the  relaxation  factor  will  mean  a  sharp  increase  in  iterations  required  for  convergence.  Although  error 
in  the  relaxation  factor  becomes  important  as  the  P3D/S  ratio  increases,  there  is  little  change  in  the  factor  itsrlf.  A 
comparison  of  Figures  B—  1 0  and  B—  1 1  reveai  that  both  give  almost  identical  relaxation  factors  for  a  particular  grid  size; 
however,  as  expected,  a  larger  number  of  iterations  is  required  for  .he  higher  P3D/S  ratio.  These  two  figures  also 
indicate  that  there  is  little  change  of  relaxation  factor  with  respect  to  the  convergence  limit  ?  ;  however,  the  number 
of  iterations  increase  as  the  convergence  limit  decreases.  Figures  B-12  and  B-13  indicate  that  the  error  decreases  very 
sharply  as  the  grid  size  decreases.  Indeed,  less  than  1%  error  may  be  achieved  with  an  8  x  8  grid.  In  developing  these 
curves,  it  was  found  that  minus  signs  appearing  in  the  output  column,  giving  maximum  change  for  each  iteration, 
indicate  the  relaxation  factor  is  too  high.  The  ideal  relaxation  factor  was  found  to  produce  only  1  or  2  minus  signs 
in  this  column.  This  feature  became  a  guide  in  arriving  at  the  best  relaxation  factor.  An  example  of  too  large  o 
relaxation  factor  may  be  found  in  the  paragraph  on  modifications,  to  follow.  Modifications  to  the  physical  conditions 
cf  the  basic  program  may  alter  the  "best"  relaxation  factor  slightly;  however,  the  curves  in  Figures  B-10  and  B-l  1  have 
been  shown  in  use  to  be  excellent  guides. 

Strain 


A  series  of  trials  determined  the  strain  relationships  ar  various  pressures  and  iength-to-width  ratios.  The  trails 
were  all  made  with  a  base  width  of  10  inches.  The  length  was  varied  irom  10  to  30  inches,  with  an  additional  series 
of  infinite  length.  The  strain  was  calculated  as  the  average  strain  of  the  centerline  parallel  to  the  narrowest  side.  The 
length  of  centerline,  used  in  this  calculation,  was  the  sum  of  the  chord  distances  between  grid  points.  Figure  3-M 
gives  the  results  of  this  series.  It  should  be  noted  that  the  curves  relating  per  cent  strain  and  d./D  ratios  for  the  case 
of  the  square  membrane,  and  the  case  of  th  1  semi-infinite  membrane,  are  very  close.  This  indicates  thot  the  strain  along 
the  centerline  of  the  membrane  is  very  close  to  that  of  the  semi-infinite  case  of  the  same  width  and  d./D  ratio;  and 
therefore,  for  all  practicable  cases,  independent  of  the  length-to-width  ratio.  This,  of  course,  would  not  be  true  for 
other  lines  across  the  membrane.  The  curves  relating  P3D/S  and  d^D  ratios  indicate  thot  tfie  center  deflection  becomes 
almost  independent  of  the  length-to-width  ratio  above  a  length-to-width  ratio  of  approximately  3: 1 . 

Modifications  to  the  Program 

The  program  for  the  membrane  has,  as  its  boundary  conditions,  fixed  edges  Ivirvg  in  the  same  plane.  Variations 
may  be  obtained  by  altering  fhes<  boundaiy  conditions.  The  boundary  may  be  given  on  arbitrary  form  by  reoding  in 
initial  values  other  than  zero.  However,  the  boundccy  must  remoin  symmetrical  about  both  axes  because  only  one- 
quarter  of  the  plate  is  considered  in  the  program.  Further  alteration  rroy  be  hod  by  allowing  the  boundaries  to  be 
supported  by  plastic  beams  which  toko  their  shape  by  responding  to  the  force  of  the  pressure  and  the  pull  of  the  supported 
membrane.  These  beams  may  also  be  placed  ocross  the  inferior  of  the  plate.  Eoges  of  the  membrone  may  be  freed  by 
reducing  the  plastic  strength  of  these  beams  to  zero.  The  pressure  need  not  be  0  canstonf  over  the  surfoce  of  '  n*  membrone. 
The  pressure  may  be  mode  a  varioble  by  introducing,  into  >he  program,  a  routine  thot  computes  the  pressure  at  eoch 
point;  or  by  reading  in  a  given  set  of  pressures  defined  at  every  point  in  ff^e  grid.  If  it  is  necessary  to  consider  the 
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entire  plate,  rather  than  just  the  quarter  plate  heretofore  considered,  the  parameters  M  and  N  following  statement  600 
should  be  changed  to  read:  M  =  2*|X  +  1 

N  =  2*|Y  +  1 

This  changes  the  indexing  so  that  the  iteration  process  covers  the  entire  plate.  By  using  combinations  of  these  conditions, 
a  wide  variety  of  boundary  and  loading  conditions  are  possible.  The  equations  have  been  isolated  from  the  iteration 
process.  They  are  listed  from  statement  155  to  the  beginning  of  the  convergence  check.  Poisson's  equation,  or  an/  othe; 
suitable  equation,  may  be  inserted  provided  the  same  parameters  are  used  to  identify  information  either  given  to  or 
required  by  other  parts  of  the  program.  Several  of  these  variations  will  be  discussed  in  the  following  paragraphs. 


Figure  B-12 


Figure  B-13 


B- 19 


Two  Parallel  Edges  Free .  Under  a  constant  pressure  looding  the  membrane  with  two  parallel  edges  free,  tikes 
on  a  cylindrical  shape  rFTat  is  independent  of  the  coordinate  perpendicular  to  the  free  edges.  For  this  condition  the 
governing  equation  is: 


In  terms  of  finite  differences: 

r^(H+i  -  +  *i.i )  p. 

The  parameters  of  this  program  care  the  same  as  those  previously  discussed,  where  applicable.  They  are  individually 
identified  in  the  comment  statements  ct  the  beginning  of  the  program.  Input  requirements  are  those  dictated  by 
READ  500  immediately  following  the  comment  statements.  Output  is  similar  for  all  of  the  programs. 
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«»  HANSEN  1R 

*  COMPILE  FORTRAN. EXECUTE  FORTRAN 

C  FUNICULAR  SHELL.  PRESSURE  CONSTANT 

C  TWO  EDGES  FREE.  CYLINDRICAL  SHAPE 

DIMENSION  2 ! 66 ) 

C  ZERO  2(1)  SAND 

1  D02I=!.66 

2  z« n=c. 

C  MAX  IY=6A 

C  . . NUMBER  OF  DIVISIONS  ALONG  DIMS 

C  ITEPL... .MAXIMUM  I T ERAT I ONS  ALLOWED 

C  CONVL... .CONVERGENCE  LIMIT. ..MAX  CHANGE  BETWEEN  ITERATIONS 

C  DIMS . VERTICAL  DIMENSION  OF  1/2  SHELL 

C  STRES... .STRESS  PER  LINEAR  DISTANCE 

C  PRESS.... PRESSURE  PER  UNIT  AREA 

C  OMEGA.... RELAXATION  FACTOR 

C  ITER..... ITERATIONS  COMPLETED 

C  Z (  I  )..... RISE  AT  POINT  I 

C  TC0N2 ... .MAX  CHANGE  BETWEEN  THIS  ITER  AND  LAST 

C  I  Pc . POINT  OF  MAX  CHANGE 

C  PS . . . PRESSURE  OVER  STRESS  RATIO 

C  STRN... c. AVERAGE  STRAIN 

C  VLAM..,,.MESH  SIZE  ALONG  DiMB 

C  INPUT 

READ  500.  !  DENT*  3  Y.  ITERL.CCNVL.DIMB.STP.ES  .PRES3. OMEGA 
500  FURMATI3 I I0f 5F10.0) 

PRINT  600 . 1  DENT 

600  FORMAT  (26HJrUMCuLAR  SHELL.  C  Yl  I  NDER/  1  8.H  I  DENT  IF  I  CAT  1  ON  ...  I  A/ ) 

PRINT  612  - 

612  FORMAT (51H  ITER  CENTER  DEFLECTION  MAX  CHANGE  IN  SHElL/) 

C  COMPUTE  ARRAY  SIZE 

N=IY+1 
Y I  -  I  Y 

YLAM=DIMB/YI 
ITEP=1 
A  T  C0N2  =  0 . 

ICSWH'-l 

1=2 

3  A -I ( i  +  l )-2.*2( l)+Z( J-l  ) 

S-Z ( 1  +  1 !-Z( 1-1  ) 

C=(4.*YLAM**2+B**2  ) **  1 .5 

D=A*8.*YLAM/C+PRES3/STK£S 

E=16.*YLAM/C 

r=Z( 1 )+D/E*OMEGA 
PS=PRES3 /STRES 

C  SET  UP  CONVERGENCE  CHECK 

ICON  1 =  T -Z ! 1 ) 

I F ( ABSF ( 1  CONI ) -CONVL)  10.10.11 

11  I CSWH=  2 

10  I F  (  ABSF (  rCONl )-AbSF i TCCN2 ) ) 12 , 12,13 
i  3  1 C0N2  =  T CON  1 
IP=  I 

12  Z ( I ) =  T 
1  =  !  +  i 

IF ( N  ~ I ) 1  CO. 14 , 3 
14  Z!  1  +  1 ) =Z  (  I  -1  ) 

GO  TO  3 

C  CHECK  CONVERGENCE 

100  GO  TO  ( 1 Ol . 102 ) » ICGwH 

C  CHECK  ITERATIONS 
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1  02  IF  (  ITEPL  —  ITER)l‘Jl  .103. IC3 
0  DATA  TC  3E  PRIM ti)  O'JT  EACH  I  T t  R A T  1  ON 

103  1  =  N 

PRINT  610, ITER. ZII)  .IP.TCON2 

610  FORMAT i I5.F18.5.13H  Z C 1 2 » 2H ) =F5. 5 ) 

C  CHECK  AVERAGE  SLOPE 

IF( .9-ZCl)/DlMBl 10*. 1GA.105 
10*  PRINT  611 

611  FORMAT (66H  AVERAGE  SLOPE  EXCEEDS  0.9  TO  1.  EQUATIONS  DIVERGE  PAST 
1 THl S  POINT) 

GO  TO  101 

C  ITEPATE  AGAIN 

105  1 TER= I TER+1 

GO  TO  * 

C  COMPUTE  STRAIN  IN  ShtLL 

Id  1=1 

ARC=0. 

1  06  APC1=S0RTFI  YLAM**2+(Z<  1  +  1 )-?.  ( I  )  )«*2  > 

APC=ARr+ARCl 
1  =  1  +  1 

IFIN-1  )99<?»  107,106 
1*7  5 TRN  =  (  A R C -D I'AU )  / •> I  •'->6  *  1 CC  . 

PRINT  6) C:  I TER.ZI ! ) » IP.TCON2 
PR  I  NT 600  *  •  DEM 


PRINT  609 .OMEGA 

6  09  FORMAT (17H  RELAX.  FACTOR. . .FI 0. 5  I 
PRINT  6.1*  P R  E  o  3 

6  Cl  FORMATI17M  PRESSURE  ........ F 1 C . 5  I 

P°INT  502  »STRE.s 

6 C2  FORMAT ( 1 SHELL  STRENGTH. .rl 0.5 ) 
PRINT  603.P5 

6'?  FORMAT  ( 1  7H  P/S . r  1  C .  5  ) 


P P I N  T  6  C  5 » v L  A  f  .D I M3 

6  05  FORMAT  C7H  L  AiVOA  Y  .  »  . FI  0.5  .2  OH  5  .  .  . ..FI  0.5} 

PR  1  NT  6 . 6 , I Y . M  R.N 

606  FORMAT  (17H  GRID  SIZE . I3.27H  STRAIN  IN  SHELL  i  ,  FI'0 . 5  ) 

PRINT  607 

607  FORMAT i ?6H.Z  0 : jT/NCES  FOR  1/2  SHELL/) 

print  cce ,  1 1  »zm  .  1  =  1  ,n» 

6  08  FORvAt  (EH  Z  (  I  2  ♦  2H  )  =Fo  •  5  *  5H  7(l?»2HJ=Fft.5»5H  Z(I2*2H)=F8.5* 

1  5  H  i  (  I  ?.2;U=F8.5»fh  7  (  I  2  ♦  2  H  )  =  F  c  »  5  ) 

C  P t  A [';  N F y  T  DATA 

CO'  TOi  ; 

o  9  r)  r  v  r 
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Hie  genera!  shape  of  the  membrane  derived  from 
this  equation  is  given  in  Figure  8-15. 

the  same  sample  problem  used  previously  was 
reused  here*  For  this  shell,  the  deflections  became  too 
large  and  the  equations  began  to  diverge.  Experience 
has  shown  that  divergence  begins  for  this  shell  at  an 
average  slope  of  about  0.9: 1 . 


Plastic  Beom  Across  Center  of  Membrane .  The 
d i fferentlai  equation  governing  the  plastic  beams  is: 


Pi 

r*  7fr*  1 


[ 


2s  Jii 

+  P’J  ■ 0 


-i*  i 


Figure  B-15.  Two  Parallel  Edges  Free 


The  second  term  is  a  load  term  giving  the  contribution  of  the  puli  of  the  membrane  to  the  shape  of  the  plastic  beam. 

The  2  in  the  numerator  crises  because  the  membrane  is  considered  symmetrical  about  the  beam.  If  the  plastic  beam  is 
placed  at  the  edge-  of  the  membrane,  without  an  adjoining  membrane  to  create  symmetry,  the  number  2  should  be  dropped. 

In  terms  of  finite  differences: 


T57? 


j 


+ 


I* 


+  u 


+  P3.V  = 


0 


Where: 


Pj  =  plastic  strength  of  the  bear*,  lbs. 

W  =  width  of  the  beam,  inches. 

The  input  requirements  are  given  by 
READ  500  and  READ  501 .  These  are  identical 
to  the  previous  requirements  except  the  second 
data  card  must  now  include  the  strength  of  the 
beam  and  the  width  of  the  beam.  Initial  values 
of  the  surface  may  be  read  in,  as  before. 

Because  the  point  of  maximum  rise  is  no  longer 
at  the  center  of  the  plate,  the  program  saves  the 
point  of  maximum  rise  by  comparing  each  point, 
as  it  is  iterated,  and  saves  the  coordinates  of  the 
point  with  the  largest  Z  value.  Strains  are  then 
computed  through  this  point  and  average  slopes 
are  considered  between  this  point  and  the  edges 
parallel  to  both  coordinate  axis.  The  point  of 
maximum  rise  is  identified  in  the  output. 

Average  strain  in  the  plastic  beam  is  also  given. 
The  large  number  of  minus  signs  in  the  maximum 
change  columns  indicate  the  relaxation  factor 
was  too  high  and  the  program  was  "overshooting." 
The  shape  of  membrane  derived  from  this 
program  is  shown  in  Figure  B-16, 


Point  of  Maximum 


Figure  B-!6.  Plastic  Beam  Across  Center  of  Membrane 


B-23 


»** 

»« 

* 


HANSEN 

PROGRAM  IS  ITERATIVE 

COMPILE  FORTRAN. EXECUTE  FORTRAN 


2R 


C  FUNICULAR  SHELL.  PRESSURE  CONSTANT 

DIMENSIONS  U»6. 66) 

C  ZERO  2 ( I . J )  oAND 

1  D02 1=1 .66 
D02 J  =  1  .66 

2  Z( I .J) =0. 


C  MAX  IX  =  MAX  IY  =  64 

C  1C0DZ.... PUNCH  COL  1C.  1  IF  Z(I.J)  GIVEN.  2  IF  NCT"GIVEN  ' 

C  IX . NUMBER  OF  DIVISION?  ALONG  DIMA 

C  IY . NUMBER  OF  DIVISIONS  ALONG  DJMli 

C  I TERL... .MAXIMUM  ITERATIONS  ALLOWED 

C  CONVL ... .CONVERGENCE  LIMIT. ..MAX  CHANGE  BETWEEN  ITERATIONS' 

C  DIMA . HORIZONTAL  DIMENSION  OF  1/4  SHELL 

C  DIMB.. . ..VERTICAL  DIMENSION  OF  1/4  SHELL 

C  SIRES. ...STRESS  PER  LINEAR  DISTANCE 

C  PRESS. ...PRESSURE  PER  UNIT  AREA 

C  PHI. .....RELAXATION  FACTOR. ..USE  1.7  FOR  AN  8  BY  8  GRID 

C  ITER . ITERATIONS  COMPLETED 

C  Z (  I  * J ) ...RIZE  AT  POINT  I  #J 

C  T C0N2 ... .MAX  CHANGE  BETWEEN  THIS  ITEP  AND  LAST 

C  IP.JP.... POINT  OF  MAX  CHANGE 

C  PS.. .....PRESSURE  OVER  STRESS  RATIO 

C  STNSX....AVE.  PERCtNT  STRAIN  OF  SHELL  THROUGH  POINT  OF  MAX  RISE 

C  AT  X  =  J,  COMPUTED  AS  IOC  TIMES  (ARC  LENGTH  -  DIMB) /DIMB 

C  STNSY....AVL.  PERCENT  STRAIN  OF  SHELL  THROUGH  POINT  OF  MAX  RISE 

C  AT  Y=  I  ,  COMPUTED  AS  ICO  TIMES  (ARC  LENGTH  - ' D I MA J /D I MA ' 

C  XLAM . MESH  SIZE  ALONG  DIMA 

C  YLAM . MESH  SIZE  ALONG  DIMS 

C  INPUT 


REAP500, JDENT,ICODZ»IX,IY,ITERL»CONVL»PHt 

500  FORMAT (5  I  10.2F1O.Q) 

READ  5C1 » D I M A  * D I M D  » S T R E S  * P R B S ? 

501  FORMAT (4F 10.G ) 

PRINT  599 ♦ i DENT 

599  F0RMAT<16H1FUNICULAK  SHELL / 1 3h  1 DLNT I F I  CAT  I  ON . . . I  3/ ) 

PRINT  600 

6CC  FORMAT  (49H  ITER  CENTER  Dt-F  LtCT  ION  MAX  CHANSc  IN  SHELL/) 
C  COMPUTE  ARRAY  .SIZE 

M=  I  X+ 1 
N  =  I Y+  1 

C  CHECK  TO  SEE  IF  INITIAL  Z(I,J)  ARE  GIVEN 

I F ( ICODZ-l ) 999,3.4 

3  READ  502, ( { Z ( I , J) ,0  =  1 ,M) , 1=1  ,N) 

5C2  FORMAT  (8Flu.O) 

4  X  I  =  I  X 
Y  1  =  I  Y 

XLAM  =  D >  M A/X I 
YLAM  =  D I  MB /Y I 
I TEP  =  I 
10  i  -2 
J  =  2 

TCON2  =0 
ICSV/H=1 
I  SV/CFf-  1 
GO  TO  100 
1  C 1  J^-.J+l 

IF  (M-J  )  1 1  ,  I  0  ,  H;0 
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12  2 ( I » J+l )  =  Z  <  l  * J-l ) 

GO  TO  100 
11  1=1+1 
J  =  2 

I  F  <  N- 1  (999*15*100 
15  I SWCH=2 

Z( I+1.J)=Z( I-l.J) 

GO  TO  100 
102  J=J+1 

IF(M-J)200, 20,21 

20  2(1  , J+l ) =Z ( 1  *  J-l ) 

2 (  1  +  1  * J ) =2 ( I-l.J) 

GO  TO  100 

21  Z ( 1+1 , J ) -I ( i-1 *J ) 

GO  TO  100 

C  COMPUTATION  PHASE 

C  COMPUTE  MEW  2(1  *  J ) 

IOC  A  =  Z! I ,J+1 )-Z (  I  *  J- 1  ) 

B  =  Z  < I  +  1,J)-ZI I-l.J) 

C  =  Z(  1  +  1  ,J)-2.*Z  (  I  »J)+ZU-1,J) 

D  =  Z ( 1 .J+l )-2.*Z ( I ,J)+Z< I ,J-1  i 

E=Z ( 1+1 , J+l ) — Z ( 1-1 . J+ll-ZII+1 . J-l >+Z ( 1-1 , J-l ) 

F  = ( 4  »  *X  L  AM*  *2 + A  *  *  2 ) +C-. 5*A*U *E  +  ( 4 . * YL AM* *  2  +  B** 2 ) *D 
G  =  PRES3 /  i  STi'cEG*2.*XLA:''i*YLAM)*<4.*XLAM**2*YLAM#*2  +  XlAM»*2*B**2 
1+YLAM**2*A**2  )  **1 .5 
H  =  2.  *(  4.*XImM**2+A**2+4.*YL A '•*>*:  5**2  ) 

T  =  Z  (  I  .  J  )  +PHI  *  (  F  +(j )  /h 
PS=PRES3/STRES 

C  SET  UP  CONVERGENCE  CHECK 

T CONI =T-Z (  I  .J  ) 

I F ( ABSF ( TCCN1 l-CONVL ) 111 ,11 1 ,107 

107  I CSWH=2 

111  I  F  (  ABSF  (  T  COM  ) -Ab^F  (  TC0N2  )  1106,106,100 

108  T COM2  =  T CON  1 
I  P=  I 

JP  =  J 

106  Z ( I , J )  =  T 

C  FIND  RETURN  S  T  A  T  t E  M  T 

109  GO  TO  (  101  .1^2)  ,  lS.-:Ch 

check  convergence 

20 C  GO  TO  (2C3.2U2) » ICjWH 
C  CHECK  ITERATIONS 

202  I  F  (  ITERL-ITLR)2'j  3, 204,204 

C  DATA  TO  BE  PRINTED  OJT  EACH  ITERATION 

204  PRINT  610,  I  TER  ,2  (N.M  *  I  P  »  JP  ,  TCON2 

610  FORMAT! 15, F16. 5, ION  7(I2,1H,I2,2H)=F9.5! 

C  CHECK  AVERAGE  sLUPE  AT  POINT  OF  MAX  RISE 

I F ( 1 ,-Z ( M  *M ) /D I Ma ) 29 , 29 , 35 

3  5  I  F  (  1  ,-Z  (  N,M  )  /  UI  i‘« A  )  29 , 29,3  0 

29  PRINT  611 

611  FORMAT (6AH  AVERAGE  oLuPt  EXCEEDS  1  TO  1,  EQUATIONS  DlVtRGL  PAST  TH 
1  I S  POINT  ) 

GO  TO  203 

C  ITERATE  AGAIN 

30  I TEP= I TER+1 
GO  TO  10 

C  COMPUTE  STRAIN  AT  X=A 

203  1=1 
J=M 
ARC  =  0 . 

31  ARC1 =SQRTF ( YLAM**2+ < L  I  1  +  1 , J ) -l (  I  ,  J )  )**2  ) 
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ARC=ARC+ARC1 
1  =  1  +  1 

I  F  ( N-  I >999.32.31 

32  STNSX= (ARC-DIMB) /DIMB*lOO. 

C  COMPUTE  STRAIN  A I  Y  =  B 

I  =  N 
J  =  1 

ARC=0. 

33  ARCl=SGRTF(XLAM**2+(Zi I  . J  +  l ) -Z ( I  . J ) ) **2 ) 

ARC=ARC+ARC 1 

J=J  +  1 

IFIM-J >999.34.33 

34  STNSV= (ARC-LMMA > /DIMA+100. 

PRINT610. ITER.Z(N.M)  .iP,JP.TCON2 
PRINT  599 . 1 DENT 

PRINT  609. PHI 

609  FORMAT ( ] 7H  RELAX.  FAC  TOR •• • F 1 C. 6 ) 

PRINT  612.N.STNSY.M.STNSX 

612  FORMAT ( 16H  oTRAIN  ALONG  1=12. 4H  IS  F3.2.18H  STRAIN  ALONG  J= I  2 . 
14H  IS  F5.2) 

PRINT  6C1.PKES3 

601  FORMAT  ( 1  ?H  PRESSURE . *10. . 

PRINT  602.STRLU 

602  FORMAT ( 1 7 H  oHtLL  ST K LNuTm. . F 1 0 . i> > 

PRINT  603. Ps 

603  FORMATI17H  P/s . FIC.b) 

PRINT  604,  XL  AM,  Ul. -'.A 

604  FORMAT  (17H  LAMDA  X.  . . F10.6 .2'0h  A . F10.5) 

PR  I  NT  6- 5  ♦  YL AM  »C I  Mb 

605  FORMAT  (  ]  7H  LA  MO  A  . F10.S.20H  . . F I  0  •  5  ) 

PRINT  606.IX.IY 

606  FORMAT  <17H  GRID  SIZE . I  ?  « 4H  SYI3) 

PRINT  607 

607  FORMAT ( 26HGZ  DISTANCES  FOR  1/4  SHLLL / ) 

PRINT  6Ce»  (  (  I  *  J  »  Z l  I  ♦  J ) »  J  = 1  » M ) » I  =  1 »  N ) 

608  FORMAT (5H  Z(I2.1H,I2»2H)=F3»S,5H  Z ( I  2 »1H , I  2  ,2H ) =F6 ,5 , 

1  5H  Z( I2»ln,l2»2ri)=F8.3.3H  Z (  I  2 . 1H , I  2  ,2H ) =F8 . 5 , 

2  SH  Z( I2,lH»I2.2H)=Fe.5»Sh  Z ( 1 2 , 1 H , I  2 ,2H ) =F8 . 5 ) 

C  READ  NEXT  L>rt T A 

GO  TO  1 
999  END 
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Pl«tle  Beams  Across  Two  Parallel  Edges.  The  same  differential  equation  governs  these  beams.  It  should  be 
n°t«d  that  the  differential  equation  considers  only  a  symmetrical  case  abou.  the  beam.  Thus,  in  this  instance,  the 
membrane  is  computed  as  though  there  were  an  infinite  set  of  such  membranns  set  end-to-end,  each  contributing  one-half 
of  the  force  carried  by  the  plastic  .earn. 

Input  requirements  are  identical  to  the  case  of  a  plastic  beam  across  the  center  of  the  membrane.  The  shape 
derived  by  this  program  is  given  in  Figure  B-17. 


Figure  B-17.  Plastic  Beams  Across  Two  Parallel  Edges 
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n  o  rm  r.  n  r*.  o  r»  o  o  o  r»  o  r*  n 


»**  HANSEN  IR 

**  PROGRAM  IS  ITERATIVE 

*  COMPILE  FORTRAN, EXECUTE  FORTRAN 


C 

C 

C 


c 

(. 

r 

c 

c 

c 


FUNICULAR  SHELL,  PRESSURE  CONSTANT 
RESTRAINING  BAND  ACROSS  X=DIMA 
DIMENSION  2(66, 66) 

ZERO  Z ( I » J  )  BAND 

1  D02 1=1,66 
D02 J  =  1  ,66 

2  Z  (  I  ,  J  )  =  0  . 


MAX  IX  =  MAX  1Y  =  64 

I CODZ ••••PUNCn  OOL  10,  1  IF  Z(I*J)  GIVEN,  2  IF  NOT  GIVEN 

IX . NUMBER  OF  DIVISIONS  ALONG  DIMA 

IY . NUMBER  OF  DIVISIONS  ALONG  Di"" 

ITERL... .MAXIMUM  ITERATIONS  ALLOWED 

CCNVL... .CONVERGENCE  LIMIT. ..VAX  CHANG;  BETWEEN  ITERATIONS 

DIMA . HORIZONTAL  DIMENSION  OF  1/4  SHFLL 

D  I  . VIRTICmL  DIMENSION  OF  1/4  SHELL 

5TRES... .STRESS  PER  LINEAR  DISTANCE 
PRE53.. . .PRESSURE  PER  UNIT  AREA 

PHI . RELAXATION  FACTOR. ..USE  1.7  FOR  AN  6  3Y  8  GRID 

ST  RNB  ••••YIlLu’  FORCE  OF  BEAM,  POUNDS 

ITER . ITtRATlONS  COMPLETED 

7.  (  I  ,  J  )  ,  .  .  R  I  L  E  A  f  POINT  I  ,  J 

T CON2 ... .MAX  Change  BETWEEN  this  1TFR  AND  LAST 
i p » jp • . • • po i n i  of  max  change 

PS . PRESSURE  OVER  STRF^S  RATIO 

s  infix  ....  avc  .  percent  strain  of  pest&aing  band  ibfamj  at  x=dima 

COMPUTED  AS  100  T I M E S  (ARC  LENGTH  -  Dliwp)/DTMB 
ST'ISX...  ,AVr.  •  PERCENT  STRAIN  0-  SE'ELL  THROUGH  POINT  OF  MAX  RISE 
AT  X  =  J,  CO'-  PUT  £D  AS  ICO  TIDE'S  (ARC  LENGTH  -  D  i  MB ) /D !MB 
S  T  N !-.  Y  .  .  .  .  A  V  l  ,  PERCENT  STRAIN  C,r  c-  HELL  THROUGH  POINT  OF  max  RISE 
AT  Y=  I  ♦  COMPUTED  A c-  10?  TI^FS  (ARC  LENGTH  -  DIMA) /DIMA 

X L A m . SIZE  ALONG  i  I*'A 

YLAM . VESrl  SIZE  ALONo  DI“ 

INPUT 


DrA0  SO'  ,  IDE  NT  ,  I  COD/.  ,  I  X  ,  I  Y  ,  I  T  f  RL  .  CONVL  ,PH  I 
SOC  F ?PM AT  (5  1  K'-»/r  1  w.C  ) 

READ  r>C  )  , D  I  T  , 0  I  m y  ,  5 T R F S  * P R r  r  ?  .ST^NB  ,W  IDTH 
5  Cl  F OR" A  T ( 6 F 1  '  .  C  ) 

PRINT  SBR.IDX'NT 

599  FORMAT  ( 4  PH  1  F  UN  l  CUL  A  K  S-.  i  L  L  L  •  r  I  >  A  1  F  I  N  G  BAND  ACROSS  X*D!MA/18H  IDE 
1 N T I F I C  A  T I  ON . . . I J / ) 


PRINT  600 

6  CO  FORMAT  (75H  i  TER  (.  L  N T  E  R  DEFLECTION  MAX  CHANGE  IN  SHELL 
;  max  CM  A  Not.  IN  uC/.V./) 

COMPUTE  array  oizf 
M  = !x+i 

N  =  I  Y  +  1 

CHECK  If'  5  F  r  l»  INI  UAL  ZU.J)  MVtN 

IF  (  lCOUZ-i 

i  1  A  0  ,  I  !  /.  I  I  ,  J  !  ,  J  ‘  1  ,  "  )  .  I  =  I  .  N  ! 

S  2  FORM/ T  I  T  F  1  .  ,  ) 

'•  x  [  r  {  X 
v  1  -  IV 

Xl.AMrDIMA/Xl 
Y  L  AM  -  0  I M  B / Y | 

I  T  L  R  -  1 
1C  1=2 
J--2 
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T  C0N2  =  U 
T  C  ON  4  =  0  , 

R  I  S  E  2  =  Q  • 

ICSWH-1 
I SWCH= 1 
GO  TO  100 
101  J*J  +  1 

I F ( IX-J) ll.10u.100 
11  1=1+1 
J--2 

1 F  ( IS!—  I  )  999, 15.100 
15  1 SWCH  =  2 

Z( 1+1.J)=Z( 1-1. J) 

GO  TO  100 
1 C 2  J=J  +  1 

I F ( IX-J)  150,21 ,21 

21  Zf 1+1 »J)=Zl 1-1 ,J) 

GO  TO  ICO 

c  computatucn  phase 
C  COMPUTE  NEW  Z ( I  ,  J  ) 

ICO  A  =  Z  (  I  ,J+1  )-Z.(  I  ,  J- 1  ) 

3  =  Z ( 1+1 ,J)-Z( I-l.J) 

C=Z( 1  +  1 , JJ-2.+Z ( I  ,J  )+Z( 1-1 , J) 

D=Z(I,J+1!-2.*Z(1.J)+Z(!,J-1) 

E=Z( 1+1 ,J  +  1 )-Z( I-i.J+1 )-Z( 1+1 » J-l )+Z ( 1-1 .J-l ) 

F= (4,*XLAM*«2+A**7 ) *C-.5*A*R*E+ ( a . * YLAM**2+B*«2 ) *D 
G  =  PRFS3/  (  STK.C  *2.  *Xl.AM*YLAV  J  *<  4.*X!_AM**2*YLAM**2  +  XLAM**2*B**2 
1  ,  YLA>/**?»A»*2  )«*!  .5 
Hr  2,*  (  4  .*XLA,V#*2  +  A**2+4.#YLAM**2  +  J**2  ) 

T=2 ( 1 »  J ) +  PH  I  * ( F  +G ) / H 
PE=PRE33/STRES 

C  SAVE  POINT  OF  f'AX  RISE 
I F ( T-R1SE2  123,23,22 

22  R I SE2  =  T 
I  05.=  I 

J  P  5.  =  J 

C  SET  UP  CONVERGENCE  CHECK 

23  T  CON  1  =  T-Z [  :  , J ) 

IF ! AUSF ( TC0N1 )-C0NVL  >11 1 .1 1  I  .10? 

127  I  _  5>WN=  2 

111  IP  i  AOSF  (  TCOM  )-AdSf  t  K0N2)  )  1  06 . 1  Of) .  1  08 

108  TcoN2=rcora 
IP- 1 

JP  =  J 

126  Z (  l  .  J  )  =T 

C  FIND  RETURN  STATEMENT 

109  GO  TO  (’01,102).! owCH 

C  COMPUTE  HE  AM 

15  2  I  *  2 
J  =  V 

155  Aft  =  Z ( I »  J- 1  I -Z l !  ,  J  ) 

CB*Z ( I* 1 , J ! -2.  +  Z ! 1 . J ) +Z  < I -1  .  J) 

3!W  (  !♦  I  .  Ji  -Z  (  I  -  1 .  J  ) 

CI**(4.»YLA,M»»2+Ut*«*2)*«1.5 

FI’*?  ,»STRE  S* AH/. n  (XI  AV»  ./.AH  »  *2  ) 

GB  *  8  •  *  6  T  /NS  «  TEAM*'  j/U; 

Gi’fWRt  Im 

Mi*  *  !  6.  •  T  RN:t *  Yi,  AV/OU 

J  (t  i  z  «  I  .  J  )  ♦  i  t  !t  /,.y  »f>H  I 

:  t  T  •  1  ON Vt  Rut  N<  .  vMtCA 

I(  ■'N',’  ■-  '  •  l  .  ^  . 


IF  tABSF  f  TC0N3 ) -CGNVL ) 151*151 » 152 
I  CSWH  =  2 

I F ( ABSF ( TC0N3 l-ABSF ( TC0N4 ) ) 153.153*154 
T  C0N4=T  C0N3 
1  PB=  1 
JP8  =  J 
Z I  I *J 1 =T8 
I  =  H  1 

IF  I N-l 1200,156,155 
Z( I+l,J)=Zt 1-1, J) 

GO  TO  155 
CHECK  CONVERGENCE 
GO  TO  (203,202 1 , ICSWH 
CHECK  ITERATIONS 
IF ( ITERL- 1  TER  1203,2 04*204 
DATA  TO  BE  PRINTED  OUT  EACH  ITERATION 
PRINT61C, ITER.ZiN.M) , I P * JP , TC0N2 , 1  PR , JPP * TC0N4 
FORMAT ( I  5  ,F 1 8 • 5  *  1CH  Z(  12,  1H,  12,  2H)=  F9.5, 

10H  Z(  12,  1H,  12,  2 H } =  F9.5) 

CHECK  AVERAGE  SLOPE  AT  POINT  OF  MAX  RISE 
1  =  IPS 
J=JPS 
XJ=J-1 

IF(1 ,-Z( I  , J 1 / (XJ*YL AM)  129,20,35 

35  X  I  =  I - 1 

IFU.-ZI  I  ,J)/(XI*XLAm;  129,29*30 

29  PRINT611 

611  FORMAT ( 64H  AVERAGE  SLOPE  EXCEEDS  1  TO  1 ,  EQUATIONS  DIVERGE  P'ST  TH 
1  IS  POINT  ) 

GO  TO  203 

C  ITERATE  AGAIN 

30  ITER^ITER+1 
GO  TO  10 

C  COMPUTE  STRAIN  IN  BEAM 

203  1=1 
J=M 

ARC  =  0. 

31  ARC1 =5QRTF ( YLAM**2+ ( Z ( 1+1 , J 1 -Z l I , J 1 1**2) 

ARC=ARC+ARC1 

1  =  1  +  1 

I  FIN- 1  1999,32.31 

32  STN«X=  (  ARC-D  I  MB  1  /  D  I  Mi?  *130. 

C  COMPUTE  A  X  I  kA  I N  IN  SHELL  AT  X-  CONSTANT 

i  =1 
J  =  JPS 
A*C=0. 

33  APC1  =  SOR  T  F  f  YLAM*  »/  +  (  Z  (  I  + 1  ,  J  1  -Z  (  I  ,  J  )  1  **2  ) 

ARC  =  ARC  + ARC  1 

1  =  1  +  1 

I  0  I  >!-  I  )  099,34,-.?. 

34  S  T N S X  =  (  ARC-!)  I  .>'i>  1  /b  I  Mb*  100. 

C  COMPUTE  MAX  SIRA  IN  IN  SHELL  AT  Y  =  CONSTANT 

I  =  I  P  S 
J  =  1 

A  P  C  =  0  . 

36  ARC1  =S()RTr  (  XlA- , .;+  (2(1  ,  J  +  l  )-Z  (  1  ,  J  )  )  *  *2  , 

ARC  =  ARC+.'iRC  I 

J=J+1 

I  r  (  M-J  1  9  9  7,  :  Uit> 

3  7  S  TNSY=  (  ARC-D  I  (Cm  )  /DIMA*  100. 

PRINTS  10  ,  I  TLR.2  IN  ,U  1  ,  I  P  ,  JP  ,  TC0N2  .  I  rp-.JI’R  ,  TC0N4 


152 

151 

154 


153 


156 

C 

200 

C 

202 

C 

204 

610 
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PRINT  599  *  I  DENT 
PRINT  609, PH  I 

609  FORMAT (17H  RELAX.  F AC  TOR . . . F 1 0 . 5 ) 

PRINT  612*IPS,STNSY, JPS.STNSX 

612  FORMAT ( 1 6H  STRAIN  ALONG  1=12, AH  IS  F5.2*18H  STRAIN  ALONG  J=12* 
1  AH  IS  F5.2  ) 

PRINT  601  ,PRES3  ,STNBX 

601  FORMAT ( 1 7H  PRESSURE . F10.5*20H  STRAIN  IN  BEAM. . .F10.S ) 

PRINT  602  sSTRES.STRNB 

602  FORMAT ( 1 7H  SHELL  STRENGTH. . FI C. 5 ,20H  BEAM  STRENGTH. .. .FI  0.5 ) 
PRINT  603, PS, WIDTH 

603  FORMAT  ( 1  7  M  P/6 . .  «F10.5»20H  BEAM  WIDTH . F10.5) 

PRINT  60A , XLAM , D I MA 

6  0A  FORMAT  { 1  7H  LA  ML)  A  X  .........  FI  0. 5 , 2  OH  A . F10.5) 

PR  I NT  635  *  YL AM  »D 1  MB 

605  FORMA  T  ( 1  7H  LAMDA  . F10.5»20H  . . F10.5) 

PRINT  606 , I  X  ,  I Y 

606  FORMAT ( 1 7H  GRID  SIZE. . . AH  6 Y  7  3  ) 

PRINT  607 

607  FORMA T ( 26HOZ  DISTANCES  FOR  1 /A  SHELL/) 

PRINT  6  C  8  »  <  (I  ,  J,Z  (  I  ,  J)  •  J=1  ,M)  ,  I  =1  ,N) 

60e  FORMAT (5H  2li2,lH,I2,£  <)-Ffi.5,5H  Z ( I  2 , 1 H , I  2 , ?H ) =F8 , 5 , 

1  5H  Z(I2,1H,I2»2H-  F8.5,5H  Z ( I  2 ♦ 1H *  I  2  ,2H ) =F8 . 5  , 

2  5  H  Zt I2,1H,12,2H*  ■  ^ . 5 , 5H  Z ( J2*1H,I2,?H)=F8,5) 

C  READ  NEXT  DATA 

GO  TO  1 
999  END 
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APPENDIX  C 


SOIL  STRUCTURE  INTERACTION 


Introduction 

The  yielding  membrane  elements,  in  addition  to  being  efficient  structural  diaphrams  in  themselves,  bring  about 
favorable  «oil-$tiuchjre  interaction  behavior.  This  Appendix  presents  background  into  the  qualitative  aspects  of  this 
behavior  and  suggests  probable  magnitudes  of  ultimate  attentuation  of  blast  overpressures  that  could  result  from  these 
interactions. 

Soil-Structure  Interaction  Forces 

Soil-structure  inreraction  forces  are  those  forces  which  act  at  the  interface  between  a  buried  structure  and  the 
surrounding  soil  medium.  These  are  generally  considered  to  be  normal  pressures  but,  to  be  sure,  shearing  forces  also  exist 
at  this  interfacial  junction. 

Normal  forces  result  from  pressures  which  act  normally  to  the  interfacial  surface  ere  generally  ?n  the  same  order 
of  magnitude  as  the  surface  overpressures  in  the  air  medium  above  the  ground  surface.  Shearing  forces  are  those  forces 
which  act  tangential  to  the  interfacial  surface  and  generally  are  in  the  same  order  of  magnitude  as  the  respective  shearing 
forces  in  the  soil  under  these  conditions.  In  general,  the  normal  forces  produce  the  grearer  effect  or  response  in  the 
structure;  because  of  this  fact,  discussion  will  be  limited  to  the  action  of  structures  under  this  interaction  component  only. 

Effective  Soil-Structure  Interaction  Pressures 

The  effective  soil -structure  interaction  pressure  is  defined  as  that  normal  pressure  distribution  which  at  a  given 
instant  will  produce  a  static  free-field  deformation  in  the  structure  equal  to  the  deformation  of  the  structure  in  the  soil 
medium  at  the  same  instant.  It  follows  then,  if  we  can  neglect  the  shearing  components  of  interaction  forces,  that  the 
moments  and  stresses  under  this  effective  pressure  will  equal  those  in  the  confined  structure  at  the  same  instant. 

Effective  soil-structure  interaction  pressures  depend  upon  the  characteristics  and  homogeneity  of  ti  e  soil,  the 
nature  of  the  loading,  and  the  stiffness  and  geometry  of  the  buried  structure  with  relation  to  the  surrounding  soil. 

Types  of  Buried  Structures 

Because  the  type  of  buried  structure  has  so  much  effect  on  the  nature  of  magnification  cr  attenuation  c'  the 
passing  overpressure,  it  is  appropriate  to  consider  these  types  in  some  detail.  The  three  basic  types,  into  which  categories 
most  buried  structures  fall,  are  the  rigid,  rigid-flexible,  and  flexible  types.  Figure  C-l  shows  examples  of  these  types 
of  structures. 

Rigid  Buried  Structures.  Rigid  buried  structures  are  those  buried  structures  which  by  definition  undergo 
negligible  deformation  upon  loading.  As  a  result  of  their  rigidity,  they  have  certain  pecuiarities  of  interaction  behavior 
which  will  be  discussed  in  more  detail  later. 

Rigid-Flexible  Buried  Structures.  A  rigid-flexible  structure  is  one  which  by  definition  exhibits  rigid 
characteristics  until  some  time  in  its  rising  loading  cycle  at  which  point  it  yields  or  flexes  in  such  a  nanner  so  as  to 
reduce  its  volume  or  alter  its  shape  considerably. 

Flexible  Buried  Structures.  Flexible  buried  structures  are  those  structures  which  by  definition  exhibit  yielding 
or  other  noticeable  structural  deforrotion  immediately  on  the  first  sign  of  an  overload.  They  continue  this  yielding  or 
reduction  in  volume  behavior  throughout  the  rising  loeding  cycle. 

Settlement  Ratios 


To  fully  realize  the  nature  of  soil-structure  interaction  phenomena,  it  U  appropriate  to  consider  three  basic 
types  of  settlement  ratios.  These  are  positive,  negative,  and  zero  ratios. 
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Figure  C-l.  Typical  Buried  Structures 
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The  Positive  Settlement  Rotio.  The  positive  settlement  ratio  by  definition  is  associated  with  the  change  in 
geometry  of  the  soil  mass  surrounding  a  buried  structure  such  that  the  structure  feels  a  vertical  load  in  excess  of  the 
resultant  of  the  dead  and  live  loads  immediately  overhead. 

The  concept  of  positive  settlement  ratio  may  be  described  visually  by  the  idec.ized  drawing  in  Figure  C-2.  It 
will  be  observed  that  the  soil  mass  surrounding  the  structure  deflects  more  than  the  vertical  column  of  soil  in  which  the 
structure  is  contained.  As  a  result  of  this  idealized  geometrical  discontinuity/  vertical  shearing  forces  are  produced  which 
add  to  the  load  that  is  normally  experienced.  The  effective  soil  structure  interaction  pressure  for  this  situation  is  then 
larger  than  that  existing  in  an  undisturbed  soil  at  this  same  depth.  To  be  sure,  this  sudden  discontinuity  does  not  usually 
exist  and  more  corbeling  action  may  be  observed/  however/  the  overall  effect  is  the  same. 

Systems  which  produce  positive  settlement  ratios  are  generally  these  which  contain  rigid  structures.  These  rigid 
structures  increase  the  overall  stiffness  of  the  vertical  column  of  soil  in  which  tltey  are  contained.  A  simplified  version 
of  this  result  is  shown  in  Figure  C-3.  If  we  assume  a  stress-strain  relationship  such  as  Hooke's  Law  to  be  valid,  the  left 
column  of  soil  will  deflect  an  amount  AL|,  where  AL]  =  pL/E.  The  column  of  soil  on  the  right,  which  contains  the 
rigid  structure,  will  deflect  an  amount  AL2  =  p(L-D)/E,  where  D  is  the  height  of  the  rigid  structure.  These  two 
displacements  differ  by  an  amount  pD/E.  The  first  is  always  abater  than  the  second,  all  other  things  being  ecjual. 

The  magnitude  of  this  difference  to  some  extent  determines  the  amount  to  which  the  pressure  reaching  the  structure  is 
increased  by  this  geometrical  action. 


Figure  C-2.  Idealized  Positive  Settlement  Ratio 
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The  Negative  Settlement  Ratio.  The  negative  settlement  ratio  by  definition  is  associated  with  the  change  in 
geometry  of  the  soil  mass  surrounding  a  buried  structure  such  that  the  structure  feels  o  verticol  load  which  is  less  than 
the  resultant  of  the  dead  and  live  loads  immediately  overhead. 

The  negative  settlement  ratio  is  shown  visually  in  Figure  C.-4.  It  will  be  observed  that  the  soil  mass  surrounding 
the  structure  deflects  less  than  the  vertical  column  of  soil  in  which  the  structure  is  contained.  As  a  result  of  this 
idealized  geometrical  discontinuity,  vertical  shearing  forces  are  produced  which  subtract  from  the  load  that  is  normally 
experienced.  The  effective  soil-structure  interaction  pressure  for  this  situation  is  then  smaller  than  that  existing  in  an 
undisturbed  soil  at  the  same  depth.  As  before,  the  sudden  discontinuity  does  not  exist  ond  in  reality  soil  arching  takes 
place  but  the  overall  effect  is  the  same. 

Systems  which  produce  negative  settlement  ratios  are  generally  those  which  contain  flexible  structures.  These 
flexible  structures  reduce  the  overall  stiffness  of  the  vertical  column  of  soil  in  which  they  are  contained.  A  simplified 
version  of  this  result  is  shown  in  Figure  05,  If  the  same  linear  stress-strain  relationship  is  assumed  as  that  previously, 
the  left  column  of  soil  will  deflect  an  amount  Alj  =  pL/E.  as  before.  The  column  of  soil  on  the  right,  which  contains 
the  flexible  structure),  will  deflect  an  amount  A  I2  ~  p(L-D)/E  +4  D.  If  AD  is  greater  than  pD/E,  then^l  >_2  will  be 
greater  than  4  Lj.  4D  will  always  be  greater  than  pD/E  if  the  structure  is  more  flexible  than  the  soil  mass  it  replaces. 
Flexible  structures  by  definition  are  not  as  stiff  as  surrounding  soil,  therefore,  they  always  produce  negative  settlement 
ratios. 


Zero  Settlement  Ratio.  A  zero  settlement  ratio  is  defined  as  being  associated  with  that  condition  when  the 
surrounding  soil  moss  and  the  soil  column  containing  the  structure  deflect  equal  amounts.  Under  such  conditions,  the 
effective  soil  structure  interaction  pressure  is  equal  to  that  existing  in  an  undisturbed  soil  medium  at  the  same  point. 
Figure  C-6  shows  such  a  condition. 


Figure  C-3.  Positive  Settlement  Ratio 
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These  systems  are  those  in  which  the  soil  and  structure  possess  equal  stiffnesses.  Certain  types  of  rigid, 
rigid-flexible,  and  flexible  structures  may  at  some  point  in  their  loading  cycle  exhibit  this  behavior.  In  general,  such 
situations  rarely  happen  throughout  the  entire  loading  cycle. 

Soil -Structure  Interactions 


Rigid  Fully-Buried  Structures.  Rigid  fully-buried  structures  generally  produce  positive  settlement  ratio 
conditions^nd  as  a  result  should  be  designed  for  pressures  in  excess  of  those  existing  at  similar  points  in  undistrubed 
soils.  By  definition,  a  rigid  structure  is  one  which  undergoes  negligible  deformation  on  loading.  According  to  the 
Air  Force  Design  Manual  (AFDM)  definition,  a  fully  buried  structure  is  one  which  is  buried  sufficiently  so  that  transient 
effects  of  shock  wave  ioadings  mcy  be  neglected.  This  arch,  if  corresponding  to  the  fully-buried  definition,  can  only 
undergo  uniform  compressive  stress  by  virtue  of  its  uniform  pressure  loadings.  The  only  bending  that  can  develop  is  due 
to  the  change  in  curvature  associated  with  the  uniform  change  in  radius  that  results  from  this  idealized  loading. 

Qualitative  aspects  of  the  behavior  of  this  structure  under  a  traveling  pressure  wave  are  shown  in  Figure  C-7. 

At  initial  contact  of  the  wavefront  with  the  structure,  non-uniform  pressures  are  developed.  These  pressures  deform 
the  cylinder  immediately  with  the  resuit  that  passive  earth  pressures  develop  on  the  sides  at  right  angles  to  the  wavefront. 
Very  rapidly,  the  situation  degenerates  or  stabilizes  into  that  shown  in  the  later  diagrams.  Once  the  pressure  is  uniform, 
it  then  starts  to  decay  somewhat  proportional  to  the  decaying  surface  wave. 

A  qualitative  load-strain  diagram  for  this  shape  is  shown  in  Figure  C-8,  Note  that  such  rigid  structures 
generally  fail  suddenly  c>n  overload. 
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Flexible  Buried  Structures.  Flexible  buried  structures  exhibit  the  opposite  behavior  to  rigid  structrues. 
Yielding  begins  almost  instantly  with  sign  of  overpressure  and  continues  until  such  has  been  relieved.  These  structures 
contain  yielding  tension  membranes  as  the  roof  and  floor.  Such  structures  will  produce  negative  settlement  ratios 
which  rapidly  attenuate  blast  overpressures. 

A  qualitative  load-strain  diagram  for  this  shape  is  shown  in  Figure  C-9.  Note  the  large  strain  which  follows 
a  low  load. 

Rigid-Flexible  Buried  Structures.  As  the  name  implies,  a  rigid-flexible  structure  exhibits  the  qualities  of 
each  type  during  its  loading  cycle.  As  might  be  expected,  positive  settlement  ratios  immediately  followed  by  negative 
ratios  may  develop.  The  qualitative  load-strain  picture  for  this  structure  may  be  seen  in  Figure  C-10. 

This  structure  is  actually  ambidextrous  in  that  it  may  exhibit  rigid,  flexible,  or  rigid-flexible  behavior 
depending  on  the  nature  of  loading,  type  of  backfill  procedure,  etc.  Generally,  however,  it  is  quite  rigid  until  either 
large  elastic  deformations  or  buckling  takes  over.  Either  of  these  latter  effects  are  those  of  a  flexible  nature. 

Figure  C-ll  shows  stable  yield  of  this  structure.  Figure  C— 12  shows  unstable  yield. 

Summary  of  Effects.  The  various  types  of  structures,  because  of  their  various  actions,  feel  different 
transmitted  pressure  waves.  These  waves,  in  their  different  forms,  may  be  seen  in  Figure  C-13.  Note  the  immediate 
advantages  of  the  rigid- flexible  and  flexible  types. 

Analysis  Features 

The  stage  has  been  set,  by  the  previous  qualitative  discussions,  for  the  statement  that  quantitative  predictions 
of  these  soil-structure  interaction  loads  are  most  difficult.  Here  we  have  a  statically  indeterminate  structural  problem 
of  the  worst  type.  Very  little  quantitative  results  of  any  kind  are  available  to  substantiate  reliable  magnitude 
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predictions.  There  is  a  particularly  intense  need  for  more  theoretical  and  experimental  data  on  the  soil-structure 
interaction  phenomenon. 

,  quantitative  predictions  are  necessary  if  hopes  of  meaningful  analyses  are  real.  Fortunately,  designs  may 

e  pr  uced  from  >vhat  limited  knowledge  we  now  have, if  we  are  not  overly  concerned  about  being  conservative. 

Design  Features 

The  ultimate  in  structural  analysis  is  to  find  an  answer,  such  as  stress  and  displacement;  given  a  structure,  its 
supporh^and  its  loads.  Obviously,  for  most  physical  systems  there  is  generally  but  one  answer.  The  analyst  hopes  to 
eir  er  ind  this  answer  exactly,  or  gain  a  close  enough  approximation  so  that  his  answer  is  acceptable.  In  short,  the 
analyst  is  a  problem  solver.  The  ultimate  in  design  is  to  create  a  given  structure,  to  resist  given  loads,  over  given 
boundaries  such  that  an  analysis  is  not  necessary  to  assure  that  this  structure  will  perform  satisfactorily.  In  short,  then, 
the  designer  is  a  problem  avoider. 

A  qualitative  understanding  of  the  general  physical  behavior  of  an  underground  structure,  as  we  have  just 
considered,  is  not  sufficient  for  analysis.  However,  such  an  understanding  ]s  sufficient  for  design.  Because  so  little 
is  known  about  the  quantitative  behavior  of  underground  structures,  as  compared  with  those  above-ground,  we 


Figure  C-6.  Zero  Settlement  Patio 
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unavoidably  find  that  our  designs  are  conservative.  This  is  not  altogether  bad,  however,  because  the  source  of  con¬ 
servatism  is  gene. ally  found  in  the  supporting  strength  offered  by  the  soil.  For  regions  in  which  blast  overpressures 
are  considered,  close-in  fallout  and  initial  radiation  will  almost  assuredly  be  such  that  quite  a  lot  of  mass  will  be 
required  for  adequate  shielding.  There  is  no  more  economical  mass  for  shielding  that  earth;  and,  therefore,  in  such 
regions,  buried  structures  make  sense  from  the  fallout  and  radiation  standpoint,  certainly  from  the  blast  resistance 
standpoint. 


UNDEFORMED  GROUND- 
SURFACE 


UNDEFORMED - 
'  STRUCTURE 


DEFORMED  GROUND 
^  SURFACE 


DEFORMED - 
STRUCTURE 


NO  2 


m 


\  ip 

;  :  3U 


r=rrrllfiJi_T 


-UNDEFORMED 
v  GROUND 
SURFACE 

j  £  >PCt) 

VUNDEFORMED 
STRUCTURE 


^-DEFORMED 
<  GROUND 
A  SURFACE 

H  i=>Pco 

—  DEFORMED  - 
STRUCTURE 


NO  5 


Figure  C-7,  A  Structure  Under  a  Traveling  Pressure  Wave 


Figure  C-8.  Qualitative  Load-Strain  Diagram  for  Rigid  Structures 
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Figure  C-9.  Qualirative  Lead-Strain  Diagram  for  Flexible  Structure! 
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Figure  C— 12-  Unstable  Yield  of  Culvert 
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Figure  C-13.  Pressure  Waves  on  Structures 
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HISTORICAL  REVIEW  OF  MEMBRANE  THEORY 

A  brief  review  of  the  development  of  large  deflection  theories  of  plates  and  the  membrane  theory  of  shells  is 
offered  in  order  to  supplement  the  previous  analyses.  The  review  is  not  intended  to  be  exhaustive  in  these  fields,  but 
rather  a  survey  pointing  out  major  contributions.  A  reasonably  complete  bibliography  supplements  the  discussion. 

In  1910,  von  Karman  extended  the  linear  theory  of  plates  by  taking  into  account  the  strain  in  the  middle  plane 
of  the  plate.  He  derived  two  non-linear  differential  equations,  as  follows: 
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In  1912,  an  exact  solution  of  von  Karman's  equations  for  a  thin, infinitely- long  rectangulai  strip  with  clamped 
or  supported  edges  was  obtained  by  Boobnov. 

In  1915,  Hencky  obtained  an  approximate  solution  o  a  laterally  loaded  circular  membrane,  or  plate  of 
negligible  flexural  rigidity,by  a  finite  difference  approach,  later  (Hencky,,  1921),  he  applied  the  same  method  to 
obtain  an  approximate  solution  of  a  laterally  loaded  rectangular  membrane.  This  same  problem  was  solved  by  F.  Foppl 
in  1922.  He  reduced  the  von  Karman  equations  to  the  following: 
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Another  approximate  solution  by  A.  Foppl  and  L.  Foppl  in  1924  made  use  of  the  Ritz  method.  They  derived 
the  following  equation  for  the  center  deflection  of  a  circular  membrane  with  clamped  edges: 
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In  1925,  Nadai  derived  an  approximate  solution  to  the  circular  membrane  problem; 
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And  in  1928  Timoshenko  used  Nadai's  approach  with  an  assumed  radial  displacement  and  derived  the  following 
equation  for  the  maximum  deflection: 
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In  1934,  S.  Way  obtained  a  power  series  of  von  Karman's  large  deflection  equations  for  the  deflection  equations, 
for  the  deflections  of  a  circular  plate.  He  carried  the  problem  out  to  a  center  deflection-to-plate  thickness  ratio  of  2/0. 
Later  (Way,  1938),  he  obtained  an  approximate  solution  for  a  clamped  rectangular  plate  using  the  Ritz  energy  method. 

In  1936,  Rudolf  Kaiser  reduced  the  general  von  Karman  equations  to  five  seuond-c.der  differential  equations; 
the  flexural  rig  dity  was  taken  as  zero  in  one  case.  These  five  equations  were  solved  numerically  by  the  use  of  finite 
difference  approximations  for  the  solution  of  a  square,  simple  supported,  plate.  The  maximum  center  deflection-to-plate 
thickness  ratio  was  2/5, 

From  a  summary  of  experimental  and  analytical  research  of  flat  plates  under  concentrated  normal  loads,  R.  G. 
Stum  and  R.  L.  Moore  (Stum  and  Moore,  1937)  concluded  that  when  the  deflection  does  not  exceed  one  half  the 
thickness  of  the  plate,  over  95  percent  of  the  load  may  be  assumed  to  be  carried  by  bending.  For  deflections  on  the 
order  of  eight  times  the  piate  thickness,  as  little  as  15  percent  of  the  load  is  carried  by  bending. 

Most  of  the  theories  of  plates  of  negligible  flexural  rigidity  I* -nit  the  maximum  deflection  to  about  two  to  four 
times  the  thickness  of  the  plate.  This  is  certainly  in  the  large  deflection  theory  range  but,  in  addition,  these  deflections 
also  cause  strains  which  are  small  enough  to  be  within  the  elastic  range. 

A  number  of  analytical  studies  have  been  conducted  to  determine  the  behavior  of  circular  membranes  (Hencky, 
1915;  Hill,  1950).  Hill  states  that,  for  the  special  case  when  the  radial  and  circumferential  strains  are  equal,  the 
strains  vary  approximately  as  the  deflection. 

Many  circular  membrane  tests  have  been  conducted  (McPherson,  et  al,  1942;  Sachs,  et  al,  1946;  Brown  and 
Sachs,  1948;  Gleyzal,  1948;  Brown  and  Thompson,  1949;  Weil  and  Newmark,  1955).  It  can  be  noted  here  that,  up  to 
the  initial  point  of  secondary  bulge  at  the  center  of  the  shell,  the  circumferential  strains  vary  approximately  as  the 
deflection  in  the  shell.  The  radial  strain  varies  from  about  one  half  the  center  strain  at  the  clamped  edge,  to  the 
maximum  at  the  center,  (at  the  center ,the  radial  and  circumferential  strains  are  equal).  These  tests  were  mainly 
interested  in  the  "instability  strains"  or  the  strain  at  the  start  of  the  secondary  bulge.  This  point  is  reached  immediately 
before  failure  occurs.  These  tests  also  indicated  that  the  circular  membranes  deflect  to  form  a  nearby  spherical  surface 
under  uniform  lateral  pressure. 

Only  two  of  tiiese  circular  membrane  tests  provided  results  which  can  be  correlated  with  the  material 
presented  here  (Gleyzal,  1948;  Weil  and  Newmark,  1955).  These  are  presented  in  Appendix  A  with  Comparison  of 
Results. 


In  an  attempt  to  determine  the  plastic  behavior  of  steel  under  a  biaxial  stress  state,  tests  have  been  run  on  thin 
tubes  (Fraenkel,  1948;  Davis  and  Parker,  1948;  Fhillips  and  Kaechele,  i956). 

in  1942,  Samuel  Levy  and  associates,  working  for  the  NACA,  conducted  a  number  of  large  deflection  tests  on 
thin  rectangular  plates  (Levy,  1942a;  Levy,  1942b;  Ramberg,  et  al,  1942).  In  one  of  these  tests  the  center  deflection 
to  plate  thickness  ratio  was  12/2.  The  results  of  this  test  are  shown  in  F'qure  D-l .  A  comparison  with  the  present 
results  is  shown  in  Appendix  A.  The  reports  deal  mainly  with  the  solution  of  von  Karman's  fundamental  equations  for 
large  deflections  by  Fourier  series.  However,  in  most  cases  the  center  deflection  of  the  plate  was  less  than  four  times 
the  piate  thickness  and  the  permanent  set  was  less  than  the  thickness  of  the  plate.  It  might  be  noted  that  they  reached 
the  conclusion  that  a  plate  with  clamped  edge;  having  c  length  to  span  ratio  of  two,  or  greater,  deflects  substantially 
the  same  as  a  plate  of  infinite  length. 

The  solution  of  von  Karman's  plate  equations  for  a  thin  membrane,  such  that  flc  oral  rigidity  is  zero,  has  also 
been  ”  ,;omplished  by  Shaw  and  Perrcne  (Shaw  and  PerTone,  1954).  They  cast  the  problem  in  terms  of  displacement 
components  u,  v,  and  vv  and  thus  obtained  three  simultaneous,  non-linear,  second-order,  partial  differential  equations. 
These  equations  were  solved  by  finite  differences  and  a  relaxation  procedure.  However,  they  dealt  with  relatively 
sn.  ill  deflections  (stresses  within  the  elastic  range). 

The  point  of  interest  here  is  that  they  solved  for  the  vertical  deflections  (w)  with  the  horizontal  displacements 
(u  and  v)  taken  to  be  zero,  and  they  also  solved  the  complete  problem  for  u,  v,  and  w.  They  noted  then  that  the  w 
displacements  ore  large  relative  to  u  and  v  and  that  the  solution  for  w  only,  holding  u  und  v  zero  everywhere,  does  not 
vary  appreciably  from  the  w  displacements  obtained  from  the  complete  solution  for  u,  v,  and  w.  That  is,  the  solution 
for  vertical  displacements  can  be  obtained  accurately  by  considering  equilibrium  in  the  vertical  direction  only. 
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In  1959,  W.  Zema  developed  the 
equations  that  must  be  solved  to  calculate  the 
stresses  in  a  membrane  shell  which  has  a 
circular  cross-section ,  at  any  point  along  its 
length  (or  width);  and  has  the  same  angle 
with  the  horizontal  at  any  point  along  its 
boundaries.  The  plan  area  may  be  quite 
irregular. 


A  series  of  tests  for  ultimate 
deflection  and  strain  has  been  conducted 
by  J.  E.  Greenspon(Greenspon,  1956). 

Only  the  results  at  failure  are  given  in  the 
comparison  with  present  results  (Appendix  A). 


In  I960,  J.  E,  Greenspon  con¬ 
sidered  the  problem  of  large  deflections 
in  a  plate  and  treated  it  as  a  membrane 
under  uniform  static  pressure  load: 


d  w  +  d  w 

<3  x2  d  y^ 


-  J  (8) 


where  P  is  the  external  lateral  pressure  and 
S  is  the  tension  per  unit  length.  For  the 
maximum  deflection  of  a  clamped 
rectangular  plate  he  obtained: 
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J 

i 

e 
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w  max 


0.164  Po2  . 


(9) 


For  the  maximum  strain  at  the  middle  of  the 
long  side  x  =  0,  y  =  b/2 


where  S  =  <r  t. 

u 


Figure  D-l.  Center  Deflection  versus  Pressure  for  f  ctangular  Plates 


The  concept  of  shell  design  by  an  inverse  procedure  has  been  suggested  by  a  number  of  investigators  (Poschl, 
1927;  Horne,  1945;  Timoshenko,  1959;  Flugge,  1960;  Harrenstien,  1961).  This  concept  is  that  of  fii  ‘ng  the  shape  to 
which  a  shell  must  adapt,  in  order  to  carry  prescribed  normal  pressure  loads  under  uniform  direct  stress.  Harrenstien 
assumed  a  uniform  compressive  membrane  force  (-S)  and  reduced  the  equilibrium  equations  of  a  membrane  shell  (Wang, 
1953)  to  one  equation 


+ 


(11) 


He  used  the  exact  expression  for  the  mean  curvature  of  a  circular  nembrare,  and  an  almost  exact  expression  for  the  mean 
curvature  of  a  rectangular  membrane,  to  solve  for  the  shape  a  membrane  must  take  to  resist  a  lateral  load  completely  by 
compressive  membrane  force. 


Harrenstien  used  four  approaches  to  solve  the  equations:  (1)  direct  Integration  to  obtain  a  closed  form  solution, 
(2)  series  solution,  (3)  numerical  approximation  by  finite  differences,  and  (4)  membrane  analogy.  Shell  structure  models 
were  built  to  conform  to  the  dimensions  predicted  by  the  solutions  of  the  equations.  The  presentation  (Harrenstien,  1961) 
illustrated  the  strength  of  membrane  shells  designed  by  this  method.  Evidence  of  this  fact  may  be  observed  by  a 
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consideration  of  the  failure  load  exhibited  by  a  0.120  in.  thick  plaster-of-paris  shell.  This  unreinforced  shell,  with  a 
diameter  of  12,0  in.,  and  a  center  loading  diameter  of  2,4  in.,  resisted  a  center  load  of  685  lbs.  before  failing. 

As  can  be  observed  from  this  review,  some  work  has  been  done  in  the  region  of  large  plastic  deflections  of 
membrane.  However,  most  of  the  work  has  been  concerned  with  circular  membranes.  Little  work  has  been  accomplished 
concerning  yielding  rectangular  membranes  supported  by  yielding  edge  beams. 
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